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This  is  the  final  technical  report  prepared  by  Calspan  Corporation 
on  a  portion  of  a  multi-phase  program  sponsored  by  the  Air  Force  Aero-Propul¬ 
sion  Laboratory,  Air  Force  Systems  Command,  Wright -Patterson  AFB,  Ohio  under 
Contract  F33615-76-C-2092,  for  the  period  1  June  1976  through  31  April  1979. 

The  work  reported  herein  was  performed  as  part  of  Phases  I  and  II  of  Project 
3066,  "Investigation  of  Rotating  Stall  and  Turbine  Heat  Transfer  in  Axial 
Flow  Turbomachinery",  with  Mr.  Marvin  A.  Stibich,  AFAPL/TBC,  as  Project  Engineer. 

Dr.  Gary  R.  Ludwig  of  the  Aerodynamic  Research  Department  of  Calspan  Corporation 
had  overall  technical  responsibility  for  the  program,  and  carried  out  the  ex¬ 


perimental  measurements.  The  theoretical  aerodynamic  and  acoustic  analyses 

v»  G  irtintlv  hv  Hre  C  am J  a  r  a j i 


TABLE  OF  CONTENTS 


SECTION  PAGE 

I  DIRECT  LIFTING -SURFACE  THEORY  FOR  THE  STEADY 

LOADING  ON  AN  ANNULAR  ROTOR  .  1 

A.  INTRODUCTION  .  1 

B.  FLOW  MODEL  AND  GOVERNING  INTEGRAL  EQUATION  .  4 

C.  SOLUTION  PROCEDURE  .  H 

D.  NUMERICAL  RESULTS  AND  DISCUSSION  .  18 

1.  Comparisons  With  2-D  Strip  Theory  .  20 

2.  Comparisons  With  Inverse  3-D  Theory  .  25 

3.  Convergence  of  Calculations  .  27 

E.  CONCLUSIONS .  28 

II  UNSTEADY  LIFTING-SURFACF  THEORY  FOR  AN  ANNULAR  ROTOR  ..  32 

A.  INTRODUCTION  .  32 

B.  FLOW  MODEL  AND  GOVERNING  EQUATIONS  .  33 

C.  SOLUTIONS  FOR  PRESSURE  MONOPOLE  AND  DIPOLE 

SINGULARITIES  .  34 

D.  PRESSURE  FIELD  FOR  ENTIRE  BLADE  ROW  .  37 

E.  UNSTEADY  INTEGRAL  EQUATION  .  39 

F.  CONCLUSIONS  .  47 


III  FURTHER  COMPARISONS  OF  APPROXIMATE  ROTOR-STATOR 

INTERACTION  THEORY  WITH  ACOUSTIC  DATA  FROM 
ANNULAR  CASCADE  . 


A.  INTRODUCTION  .  48 

B.  ROTOR-STATOR  INTERACTION  NOISE  MEASUREMENTS  .  49 

1.  Rotating  Annular  Cascade  Facility  .  49 

2.  Acoustic  Measurements  .  50 


v 


SECTION 


PAGE 


C.  COMPARISONS  WITH  THEORY  .  59 

1.  Theoretical  Model  and  Input  Data  .  59 

2.  Discussion  .  60 

D.  CONCLUSIONS  .  65 

IV  SUNWARY  .  67 

APPENDIX  A.  EQUIVALENCE  OF  PRESSURE  DIPOLE  AND 
HORSESHOE  VORTEX  REPRESENTATIONS 
OF  THE  STEADY  LOADING  ON  AN 

ANNULAR  BLADE  ROW  .  71 

APPENDIX  B.  EVALUATION  OF  RADIAL  INTEGRALS  .  79 

APPENDIX  C.  EVALUATION  OF  AXIAL  INTEGRALS  .  81 

REFERENCES  .  12g 


vi 


LIST  OF  ILLUSTRATIONS 


FIGURE  PAGE 

1  Blade  Row  and  Duct  Geometry  .  87 

2  Camber  Line  Geometry  in  Blade-Fixed  Coordinates  .  88 

3  Comparison  of  Present  3-D  Theory  and 

2-D  Strip  Theory  for  Cases  1-3  . . .  89 

4  Comparison  of  Chordwise  Load  Distributions  at 

Mid-Annulus  .  90 

5  Comparison  of  Chordwise  Load  Distribution  at 
Mid-Annulus  for  Case  3  With  Isolated  Flat 

Plate  Distribution  .  91 

6  Comparison  of  Present  3-D  Theory  and 

2-D  Strip  Theory  for  Case  4  .  92 

Comparison  of  Chordwise  Load  Distributions  at 

Hub,  Mid-Annulus,  and  Tip  Radii  for  Case  4  . .  93 

8  Comparison  of  Present  Theory  with  Design  Values  of 

Okurounmu  and  McCune  for  Case  5  .  94 

9  Comparison  of  Chordwise  Load  Distributions  at 

Mid-Annulus  for  Case  5  .  95 

10  Comparison  of  Present  Theory  With  Design  Values  of 

Okurounmu  and  McCune  for  Case  6  .  96 

11  Convergence  of  Calculations  vs.  Number  of 

Azimuthal  Modes  for  Cases  1  and  4  .  97 

12  Convergence  of  Calculations  vs.  Number  of 

Azimuthal  Modes  for  Case  5  . 98 

13  Convergence  of  Calculations  vs.  Number  of 

Radial  Modes  for  Cases  1,  4  and  5  .  99 

14  Test  Configuration  for  Acoustic  Studies  of 

Rotor-Stator  Interaction  .  100 

15  Acoustic  Test  Instrumentation  .  101 

16  Oscilloscope  Records  of  Noise  Signal  from 

Rotor-Stator  Interaction  .  102 


vii 


LIST  OF  ILLUSTRATIONS  (Cont'd.) 


FIGURE  PAGE 

17  Uncorrected  1/3  Octave  Spectrum  of  Outer  Wall 
Sound  Pressure  Level  from  Rotor-Stator  Interaction 

(a)  Rotor  RPM  =  1140  .  103 

(b)  Rotor  RPM  =  1150  .  104 

(c)  Rotor  RPM  =  1160.S  .  105 

18  Uncorrected  1/3  Octave  Spectrum  of  Outer  Wall 
Sound  Pressure  Level  from  Rotor-Stator  Interaction 

(a)  Rotor  RPM  =  1695.5  .  106 

(b)  Rotor  RPM  =  1700.5  .  107 

(c)  Rotor  RPM  =  1710.5  .  108 

(d)  Rotor  RPM  =  1716  .  109 


19 


20 


21 


24 


25 


26 


Background  Noise  Correction  Applied  to 

Measured  Sound  Pressure  Levels  . . .  110 

Measured  Frequency  Response  of  Multifilters  in 

Real  Time  Analyzer  .  Ill 

Frequency  Response  of  Condenser  Microphone  .  112 

RMS  Wall  Pressure  at  Blade  Passage  Frequency, 

Stator  Stagger  Angle  £  c.=  37.2  deg  .  113 

JIM 

RMS  Wall  Pressure  at  Second  Harmonic  of 
Blade  Passage  Frequency,  Stator  Stagger  Angle, 
h  SM  =  37.2  deg  .  114 

RMS  Wall  Pressure  at  Fourth  Harmonic  of 
Blade  Passage  Frequency,  Stator  Stagger  Angle, 

§SM  3  37 . 2  deg  .  115 

RMS  Wall  Pressure  at  Fifth  Harmonic  of 
Blade  Passage  Frequency,  Stator  Stagger  Angle, 

<5  SM  =  3  '  “  .  116 

RMS  Wall  Pressure  at  Blade  Passage  Frequency, 

Stator  Stagger  Angle  *  28.2  deg  .  117 

RMS  Wall  Pressure  at  Second  and  Third  Harmonics  of 
Blade  Passage  Frequency,  Stator  Stagger  Angle, 

SsM  ’  -8-2 


118 


FIGURE 

:s 

29 


30 


31 


LIST  OF  ILLUSTRATIONS  (Cont'd.) 


RMS  Wall  Pressure  at  Fourth  Harmonic  of 
31ade  Passage  Frequency,  Stator  Stagger  Angle, 
SSM  1  23.:  deg  . 

RMS  Wall  Pressure  at  Fifth  Harmonic  of 
Blade  Passage  Frequency,  Stator  Stagger  Angle, 
SSM  31  23.:  deg  . 

Comparison  of  Theoretically  Predicted  SPL 
With  Experiment;  Stator  Stagger  Angle  =  23.2° 

(a)  Blade  Passage  Frequency  . 

(b)  2nd  and  3rd  Harmonics  . 

(c)  4th  Harmonic  . 

id)  Sth  Harmonic  . 

Comparison  of  Theoretically  Predicted  SPL 
With  Experiment;  Stator  Stagger  Angle  =  3'. 2° 

la)  Blade  Passage  Frequency  . 

(b)  2nd  Harmonic  . 

ic)  4th  Harmonic  . 

id)  Sth  Harmonic  . 


LIST  OF  TABLES 


PAGE 


Input  Parameters  for  Cases  1-6  .  19 

Sample  of  Acoustic  Data  Reduction  for  Spectrum  Kith 

Multiple  Harmonics  of  Blade  Passage  Frequency  .  56 

Sample  of  Acoustic  Data  Reduction  for  Spectrum  With 

Two  Adjacent  Third-Octave  Bands  Responding  Nearly 

Equally  to  the  Same  Harmonic  of  Blade  Passage  Frequency..  57 


x 


NOMENCLATURE 


CL  sound  speed 
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radial  eigenvalue 
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Leer)  sectional  lift  per  unit  span 
L&  tangential  blade  spacing,  2TCr/ & 
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NOMENCLATURE  (Cont'd.) 


Mc<n  sectional  pitching  moment  about  mid-chord,  positive  clockwise 

term  at  which  i,  j,  k,  and  m  series  are  truncated,  respectively 

NK,  Nn 

NL, NN  number  of  chordwise  and  spanwise  collocation  points,  respectively 
/V|ct  function  describing  the  blade  camber  line 


p  perturbation  pressure;  also  used  in  Section  II  to  denote  integral 
shift  in  azimuthal  mode  number  caused  by  inter-blade  phase  shift, 
Eqn .  (49b) 

-a  1  2 

p  nondimensional  perturbation  pressure,  p/  p_co  U 

cylindrical  polar  coordinates  (Fig.  1) 


^mBk  orthonormal  radial  eigenfunctions 


S, n  local  streamwise  and  normal  coordinates  (Fig.  2) 

A 

S  nondimensional  streamwise  coordinate,  Eqn.  (10) 


sign  function,  ±  f  for  X  k  0 

-i 

axial  sine  integrals  defined  in  Eqns .  (14d,  h) 

U  uniform  axial  velocity 

Up  relative  free 

qr  perturbation  velocity 

Vrn&k  defined  in  Eqn.  (2c) 

X  nondimensional  axial  coordinate,  Eqn.  (5b) 

OC  angle  of  attack 

(  7  +  M  2  )  *and  (  1  -  )  /a  ,  respectively 

r  sectional  circulation,  positive  counter-clockwise  in  Fig.  2 
<5t%)  Dirac  delta  function 

Ag , AD  turning  angles  in  blade  and  duct-fixed  coordinates,  respectively; 
Eqns.  (29)  and  (30) 

5  helical  coordinate,  &  -  ~p  Z 


■  tor*  v® 


-stream  velocity,  U  '\fl  +  \~Tj' 


NOMENCLATURE  (Cont'd.) 


% 

\ v,0k 
M-fl 

P 

<T 

<P 

<*r 

V 

60 

n 

H 

T 

o 

t  oo 


c«./2rT 

defined  in  Bqn.  (2b) 

defined  in  Eqn.  (2a)  for  steady  analysis  of  Section  I;  defined  in 
Eqn.  (44)  for  unsteady  analysis  of  Section  II. 

inter-blade  phase  shift,  Eqn.  (49a) 

density 

nondimensional  radial  coordinate,  r /  rr 
cm'1  % 

U>rT  /  U 

undisturbed  flow  angle,  tan 

rotor  angular  velocity  in  steady  analysis  of  Section  I;  harmonic 
excitation  frequency  in  unsteady  analysis  of  Section  II. 

rotor  angular  velocity  in  unsteady  analysis  of  Section  II. 


evaluated  at  hub 
evaluated  at  tip 


source  coordinate  on  surface,  e.g.,  X0  ,  C~0  %  also  used  in 
Section  II  to  denote  conditions  at  undisturbed  state. 


conditions  at  2  =  ±  <* 


(  ) 


Superscripts 


dimensionless  quantity 
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SECTION  I 


DIRECT  LIFTING-SURFACE  THEORY  FOR  THE 
STEADY  LOADING  ON  AN  ANNULAR  ROTOR 

A.  INTRODUCTION 

Increased  emphasis  on  reduction  of  the  size,  weight,  and  noise  output 
of  axial  flow  turbomachinery  demands  improved  understanding  of  the  flow 
through  high-speed  fan  and  compressor  blade  rows.  The  task  of  calculating 
the  fully  nonlinear,  three-dimensional,  viscous  flow  through  a  blade  row  is  a 
formidable  one  indeed.  Consequently,  some  approximations  are  required  in  order 
to  obtain  a  tractable  model,  the  most  familiar  being  the  idealization  of  inviscid 
flow  through  a  two-dimensional  cascade.  However,  as  more  detailed  questions 
are  asked  about  modern  blade  row  performance,  the  essentially  three-dimensional 
character  of  the  flow  takes  on  increased  importance. 

This  section  is  a  report  on  our  development  of  a  three-dimensional, 
direct  lifting-surface  solution  for  the  linearized,  inviscid,  steady  flow  through 
a  rotor.  Specifically,  it  presents  a  collocation  procedure  for  inverting  the 
linear  integral  equation  which  relates  the  local  loading  on  the  blades  to  a 
prescribed  camber  line;  this  integral  equation  had  been  derived  as  part  of 
our  effort  under  a  previous  contract. *  Although  such  an  analysis  is  in  principle 
limited  to  low  pressure  ratios,  it  holds  the  promise  of  increasing  our  physical 
understanding  of  the  rotor  flowfield  at  a  cost  well  below  that  of  numerically 
solving  the  full  nonlinear  equations.  The  present  approach  also  contains  impor¬ 
tant  features  not  present  in  the  two-dimensional  approximation.  For  example, 
it  does  include  disturbances  induced  by  the  trailing  vortex  wakes  which  result 
from  variations  in  the  blade  circulation,  and  the  results  reported  below  show 
this  to  be  an  important  effect.  The  analysis  and  techniques  developed  here 
should  also  prove  useful  in  the  study  of  three-dimensional,  unsteady  loading, 
treated  in  the  next  section. 


1 


The  study  of  linearized,  compressible  three-dimensional  How  through 
an  annular  blade  row  began  with  the  early  work  of  McCune,  who  initials  con- 

2-4 

sidered  only  the  disturbance  generated  by  the  blades’  thickness  distribution. 
Okurounmu  and  McCune'1  extended  the  analysis  to  treat  the  loading  problem  wherein 
the  circulation,  PiO  ,  is  specified.  Initially,  they  modeled  each  blade  as 
a  radially-oriented  line  vortex  of  variable  spanwise  strength,  and  were  able 
to  obtain  useful  expressions  for  the  overall  static  and  total  pressure  rise,  and 
change  in  axial  velocity.  They  later  generalized  this  work  to  include  finite 
chord  effects  in  a  lifting-surface  type  analysis,  with  the  blades  modeled  as  a 
distribution  of  these  radially-oriented  vortices. b  Again,  only  the  design 
problem  was  considered,  i.e.,  computing  the  camber  line  required  to  produce  a 
prescribed  loading;  this  can  be  reduced  to  an  integration  over  the  known  vortex 
distribution. 

Just  as  in  isolated  airfoil  theory,  the  solution  to  the  off-design 

problem,  in  which  one  seeks  the  loading  produced  by  blades  with  a  prescribed 

camber  line,  is  more  difficult  because  it  requires  the  inversion  of  an  integral 

£ 

equation.  McCune  and  Dharwadkar  considered  this  problem  in  the  context  of 

a  Prandtl  lifting-line  approximation,  using  the  vortex  type  of  analysis  presented 

in  Ref.  5.  This  work  represented  the  extension  to  subsonic  flow  of  Falcao’s 

earlier  investigation  of  the  incompressible  case.''  Nlamba^’^  was  the  first  to 

treat  the  off-design,  or  direct,  problem  using  a  lifting-surface  type  of 

analysis.  He  used  pressure  dipole  singularities  to  represent  the  blade  loading, 

in  contrast  to  the  vortex  representation  of  McCune,  et  al.  '  Sal  nun  later 

applied  a  pressure  dipole  representation  to  the  prediction  of  the  unsteady  blade 
1  ’  15 

loading.  ’  Although  not  concerned  with  the  steady  problem  per  se,  Salaun 
briefly  discussed  it  in  Appendix  P  of  Ref.  15;  however,  no  numerical  results 
were  obtained  for  this  case. 

Thus,  the  present  work,  which  also  employs  a  pressure  dipole  approach, 
most  closely  parallels  that  of  Namba;  in  fact,  the  basic  flow  model  and  the 
assumptions  made  are  almost  identical.  However,  we  found  that  our  form  of  the 
resulting  integral  equation  differed  from  Mamba's  in  several  significant 
respects.  Accordingly,  a  major  portion  of  our  effort  under  the  previous  contract 


was  spent  in  trying  to  reconcile  these  differences,  or  barring  that,  to  decide 
which  formulation  was  correct.  We  undertook  a  detailed  review  of  all  the  steps 
leading  to  the  formulation  of  the  integral  equation.  In  particular,  it  was 
verified  that  our  pressure  monopole  and  dipole  singularity  solutions  do,  in 
fact,  satisfy  the  governing  partial  differential  equation,  and  that  our  expressions 
for  the  pressure  and  velocity  fields  associated  with  the  entire  blade  row 
exhibit  the  appropriate  discontinuities  across  the  blade  surfaces  and  trailing 
vortex  wakes.  It  was  also  shown  that  the  pressure  and  velocity  fields  satisfy 
global  mass  and  momentum  balances  in  the  duct.  On  this  basis,  it  was  concluded 
that  our  formulation  of  the  integral  equation  is  the  correct  one. 

This  review  of  the  linearized  analysis  was  reported  in  detail  in  Ref.  1. 

During  the  current  contract,  we  have  also  found  that  our  formulation  of  the 

problem  agrees  with  that  given  by  Salaun.  In  addition,  Salaun's  ingenious 

manipulation  of  the  infinite  series  which  arise  provided  us  with  the  link 

necessary  to  demonstrate  the  analytical  equivalence  between  the  horseshoe 
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vortex  representation  of  the  blades  used  by  McCune  et  al.,  and  the  pressure 
dipole  representation  used  here.  This  equivalence  is  discussed  in  Appendix  A, 
ind  further  supports  the  validity  of  the  present  integral  equation. 

The  major  effort  under  the  current  contract  has  been  spent  in  developing 
an  efficient  means  of  inverting  the  integral  equation  governing  the  direct 
lifting-surface  problem.  As  will  be  seen  below,  the  method  of  solution  is  the 
second  feature  distinguishing  the  present  investigation  from  that  of  N'amba.* 

The  procedure  adopted  here  is  felt  to  be  a  more  direct  extension  of  the  kernel- 
function  methods  frequently  used  in  isolated  airfoil  theory.  Moreover,  the  form 
assumed  for  the  unknown  loading  allows  analytical  expressions  to  be  derived  for 
the  required  spanwise  and  chordwise  integrations,  thus  affording  greater  accuracy 
and  improved  efficiency  in  terms  of  computer  time. 

Sub-section  B  briefly  reviews  the  basic  flow  model  and  introduces  the 
governing  integral  equation.  Sub-section  C  describes  the  solution  procedure, 
ind  numerical  results  and  comparisons  with  other  theories  are  presented  and 
discussed  in  Sub-section  D. 


3 


B. 


FLOW  MODEL  AND  GOVERNING  INTEGRAL  EQUATION 


The  model  assumes  the  blade  row  is  housed  in  an  infinitely  long,  hard- 
walled  annular  duct  of  constant  hub/tip  ratio,  h  ,  containing  a  uniform  sub¬ 
sonic  axial  flow  at  Mach  number  M  ,  as  shown  in  Fig.  1.  The  flow  is  assumed 
to  be  inviscid  and  steady  in  blade-fixed  coordinates.  In  this  reference  frame, 
the  inflow  has  a  relative  velocity  =  [U2+(oj^)*J  2  ,  and  follows  the 
helical  undisturbed  stream  surfaces  defined  by  z,  *  9  -  2  =  constant. 

The  blades  are  assumed  to  produce  only  small  perturbations  about  this  undis¬ 
turbed  flow,  allowing  the  linearised  boundary  conditions  to  be  applied  along 
these  surfaces.  The  inviscid  and  small  perturbation  assumptions  further  imply 
that  the  disturbance  flowfield  is  irrotational ,  and  hence  isentropic.  We 
restrict  ourselves  here  to  considering  only  the  case  where  the  undisturbed 
relative  Mach  number  at  the  tip  is  subsonic.  This  ensures  that  all  disturbances 
will  damp  out  far  upstream  of  the  rotor,  and  that  far  downstream  only  those  aue 
to  the  trailing  vortices  will  remain. 

Since  the  problem  has  been  linearised,  the  disturbance  fields 
produced  by  blade  thickness  and  loading  may  be  determined  separately  and  then 
superimposed.  However,  this  separation  is  not  as  straightforward  as  it  is  for 
a  planar,  isolated  airfoil.  In  the  thickness  part  of  the  rotor  problem,  the  blades 
must  be  cambered  to  achieve  the  condition  of  tero  loading.  This  sero- loading 
camber  line,  or  "thickness- induced  camber",  can  be  calculated  as  an  integral  over 
the  known  source  distribution  once  the  blade  thickness  and  the  operating  condi¬ 
tions  are  prescribed.  This  point  is  discussed  at  more  length  in  Refs.  1  and  14; 
the  latter  contains  sample  calculations  as  well.  In  the  ideal  iced  case  of 
infinitely  thin  blades,  the  present  loading  analysis  can  be  applied  directly 
without  consideration  of  this  effect.  However,  for  blades  of  finite  thickness, 
the  analysis  assumes  that  the  thickness- induced  camber  has  been  calculated  and 
that  the  geometric  camber  is  measured  from  this  zero- load ing  camber  line. 


As  mentioned  earlier,  the  derivation  of  the  governing  integral  equation 
appears  in  a  previous  report,1  and  so  will  not  be  repeated  here  except  for  a 
brief  outline  to  introduce  some  concepts.  The  first  step  is  in  the  manipulation 
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of  the  linearized  flow  equations  into  a  single  equation  for  the  perturbation 
pressure,  -p  .  A  singular  solution  to  this  equation,  or  Green's  function,  is 
then  found  which  physically  represents  an  idealized  point  force/volume  acting 
on  the  fluid.  This  will  be  referred  to  as  the  pressure  dipole  solution. 


To  obtain  the  pressure  field  for  the  entire  blade  row,  these  singulari¬ 
ties  are  distributed  along  the  undisturbed  stream  surfaces  c,  -  "5*“’  ft  =0, 1  •  g-  1  , 
with  the  dipole's  strength  equal  to  the  local  loading,  A-p  cra  ,  s0 )  .  On  a 
blade  surface,  is  uniquely  related  to  2a  ,  and  so  we  set  A  p  ( <~a  ,  s.d, ) ) 

-  Ap(f"0  ,  2a  )  ,  and  transform  the  chordwise  integration  to  one  over  the  axial 
coordinate.  The  resulting  expression,  Eqn.  (134)  in  Ref.  1,  is: 
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and  the  remaining  symbols  are  defined  in  the  Nomenclature.  It  has  been  assumed 
tor  convenience  in  Eqn.  (1)  that  the  axial  projection  of  the  blade  chord,  C ^ 
in  Fig.  2,  is  independent  of  radius,  though  this  is  not  essential  to  the  analysis. 
Eqn.  (1)  has  been  shown  to  yield  a  pressure  field  which  is  continuous  everywhere 
except  across  the  blades,  and  satisfies  the  same  hard-wall  boundary  condition 
Imposed  on  the  pressure  dipole  solution, 
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Eqn.  (.3)  is  equivalent  to  requiring  that  the  radial  velocity  vanish  at  the  hub 
and  tip. 


The  above  result  is  a  superposition  of  duct  modes  which  are  sinusoidal 
in  d  with  a  fundamental  period  -g—  .  have  a  complex  exponential  dependence 
on  2  ,  and  whose  radial  variation  is  described  by  the  radial  eigenfunctions 

and  their  associated  eigenvalues,  K^6h  .  The  denote  linear  combi¬ 

nations  of  Bessel  functions  of  the  first  and  second  kind  which  satisfy  liqn.  (.3). 

Their  argument  is  actually  K\^gkr  .  which  is  written  here  simply  as  cr  for 

"'-4 

the  sake  of  brevity.  These  functions  are  described  at  length  by  McCune'  and 

,  .  .  15 

Tvler  and  sot r in. 


knowing  -p  ,  the  normal  component  of  the  disturbance  velocity,  , 

(.see  Fig.  2)  can  be  obtained  by  integrating  the  corresponding  momentum  equation 
from  upstream  infinity  along  the  undisturbed  helical  streamline  through 
(r,  9  ,  j2)  .  This  result,  obtained  by  combining  Eqns.  (.145)  and  (164)  of  Ref. 


mgkv  »  »*« 


a  4>(r. 


Note  that  the  terms  in  the  first  four  lines  of  Eqn.  (.4)  are  identically 
tero  upstream  of  the  blade  row.  The  first  two  lines  represent  contributions 
from  the  no»c>  ,  k  •  Q  ,  or  plane  wave  mode.  The  next  term,  which  involves 
radial  derivatives  of  the  loading,  represents  the  influence  of  the  trailing 
vortex  pattern.  For  fixed  C,  ,  this  term  is  independent  of  £  downstream  of 
the  blade  row,  which  reflects  the  fact  that  the  vortices  are  convected  downstream 
without  attenuation  along  the  undisturbed  streamlines.  This  will  be  referred  to 
below  as  the  wake  term. 

In  the  last  group  of  terms,  since  we  have  assumed  <  i  everywhere, 

one  can  easily  show  that  from  Eqn.  (2a)  is  always  real  and  positive. 

itence,  these  terms  rapidly  attenuate  as  one  moves  away  from  the  blade  row,  and 

will  be  referred  to  below  as  the  exponentially  decaying  terms.  When  the  tip 

relative  Mach  number  becomes  supersonic,  certain  modes  will  propagate  undamped 
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down  the  duct.  This  was  first  noted  by  McCune,  who  was  thus  able  to  relate 
the  resulting  acoustic  radiation  directly  to  the  blade  thickness  distribution. 
Such  modes  are  identical  in  form  to  the  propagating  duct  acoustic  modes  studied 
by  Tyler  and  Sofrin,^  which  have  since  formed  the  basis  for  numerous  investiga¬ 
tions  of  turbomachinery  noise. 

The  above  expression  for  vn  was  shown  in  Ref.  1  to  be  continuous 
across  the  blade  surfaces,  as  is  required  physically.  As  part  of  this  demon¬ 
stration,  it  was  found  useful  to  integrate  the  exponentially  decaying  terms  by 
parts  with  respect  to  £0  .  At  that  time,  it  was  anticipated  that  this  altered 

form  for  ,  Eqn.  (172)  of  Ref.  1,  would  be  used  as  the  basis  for  the  desired 

integral  equation.  This  approach  was  later  abandoned  for  several  reasons,  the 

d&P 

principal  one  being  the  added  complication  that  appears  in  the  integrand 

rather  than  <A-p  itself,  and  the  former  can  be  expected  to  exhibit  singularities 
at  both  the  leading  and  trailing  edges. 

To  convert  Eqn.  (4)  into  the  desired  integral  equation,  the  right-hand 
side  is  specialized  to  the  reference  blade  surface,  ^  =  O  .  It  is  also  conven¬ 
ient  for  machine  calculations  to  rationalize  the  resulting  expression  so  that 
only  real  arithmetic  is  involved.  From  here  on,  we  will  work  with  the  following 
dimensionless  coordinates, 
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One  can  express  the  aspect  ratio  and  a  solidity  parameter  based  on  the  constant 
axial  chord  projection,  C*,  ,  in  terms  of  these  dimensionless  groups: 
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The  loading  is  nondimensional i;ed  by  the  dynamic  pressure  associated  with  the 
uniform  axial  velocity, 
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The  left-hand  side  of  Eqn.  (4)  is  related  to  the  blade  geometry  as 
follows.  Let  denote  the  function  describing  the  blade  camber  line 

as  a  function  of  streamwise  and  radial  coordinates  (see  Fig.  2),  after  accounting 
for  any  thickness-induced  camber  as  noted  above.  Then  the  blade  boundary  condi¬ 
tion  requires  that 

vjs-o)  -  *  J°  ^ 

where  we  have  nondimensionalited  ^C(_  and  S  by  the  local  semi-chord: 

£  s  S-c/z 

Tu  2  m  ~  (ld) 

On  the  reference  blade  2  -  S  Jr  ,  where  Iff  =  tan  *  (  <pr<r)  is  the  local  un¬ 
disturbed  flow  angle  (Fig.  2) ,  and  =  C  y  .  It  then  follows  from  the 

definitions  in  Eqns.  (5)  and  (10)  that  S-X. 

With  the  above  definitions,  the  integral  equation  relating  the  blade 
camber  line  to  the  unknown  loading  can  be  written  as: 
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To  facilitate  the  solution  procedure  described  next,  like  powers  ot  fi  have 
been  grouped  together  in  the  exponentially  decaying  terms.  The  coefficients 
of  (To  '  in  these  terms  are  given  by: 
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As  noted  in  the  Introduction,  and  discussed  more  extensively  in  Ref.  1, 
tins  integral  equation  does  not  agree  with  that  derived  by  Namba.  The  differences 
can  be  traced  back  to  the  solution  for  the  pressure  field,  Eqn.  (.1).  Samba's 
expression  for  -p  (Ref.  11,  Eqn.  (15)),  quotes  a  different  form  for  the 
rn:o  ,  k  -o  term,  and  also  contains  a  "scale  factor",  in  his 

terminology,  which  does  not  appear  in  our  result. 

Before  proceeding  to  solve  Eqn.  (11),  the  singularities  which  can  arise 
in  the  integrand  warrant  some  discussion.  In  the  form  in  which  the  wake  terms 
first  appeared,  Eqn.  (161)  of  Ref.  1,  the  rr>  and  K  series  could  be  shown  to 
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diverge  when  ■  T  and  x,  <  X  .  An  analogous  singularitv  occurs  in  isolated 

airfoil  theory.  Ref.  lo,  Eqn.  v'-CSa);  both  are  produced  b>  the  trailing  vorticity. 
\s  pointed  out  in  Ref.  lt>.  an  integration  by  parts  along  the  span  will  reduce 
the  order  of  the  singularity  so  that  it  can  be  evaluated  using  a  Cauchy  princi¬ 
pal  value.  In  Ref.  1.  this  process  was  carried  one  step  further  b>  integrating 
by  parts  twice  along  the  span,  and  making  use  of  Eqn.  lo)  at  the  limits.  This 
leads  to  the  form  shown  above,  which  involves  only  uniforml>  convergent  series. 

The  price  one  pays  is  that  now  the  integrand  involves  radial  derivatives  of  the 
unknown  dp  ;  however,  this  is  easily  handled  by  the  present  solution  procedure. 


In  the  last  group  of  terms, 
assured  by  the  exponential  decay  wit 
absent  at  x^  *  x  .  where  the  terms 
vanish  identically,  but  the 
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terms  form  divergent  series.  However,  the 
is  through  the  a  -  Xa )  ,  so  that  the 

preted  in  the  Cauchy  principal  value  sense, 
from  the  bound  vorticity  at  (  X0  ,  S~,  ) 
analytical  expression  for  the  integrals  can 
pretation  is  automatically  satisfied.  This 
and  in  Appendix  C. 


onlv  dependence  on  x,  that  remains 
axial  integrations  can  be  inter- 
Physicall'.  this  singularit'  arises 
As  pointed  out  in  Ref.  lo.  if  an 
be  found,  the  principal  value  mter- 
is  the  case  here,  as  shown  below 


C.  SOLUTION  PROCEDURE 

In  order  to  solve  the  integral  equation,  it  vs  assumed  that  the  unknown 
loading  can  be  expanded  in  a  finite  series  of  suitable  chosen  chordwiso  and 
spanwise  loading  functions.  Specit'ical ly ,  we  set 
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where 


COO 


and  the  ,  .  *  t,2  .  ..NI,  -j  -  1,2  NJ,  are  constants  to  be  determined 

by  the  solution.  The  chordwise  functions  are  the  same  as  those  often  used  in 
subsonic  airfoil  theory.  They  automatically  satisf>  the  kuttu  condition  at  the 
trailing  edge,  xt  -  *  1  ,  and  the  .  *  1  term  exhibits  the  appropriate  square 

root  singularity  at  the  leading  edge.  -1  .  The  pol\-nomial  variation  in 

<T#  allows  simple  treatment  of  the  boundary  condition  in  Eqn.  (o',  as  well  as 
the  radial  derivatives  in  the  wake  term,  as  will  be  seen  shortly. 


Substitution  of  Eqn.  (Id)  into  Eqn.  (11a)  leads  to  radial  integrations 
of  the  form: 
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the  axial  integrations  which  result  are  more  complex;  for  <.  -  •  they  are: 


'  /  -  x. 


V  >  *  X. 


t  a  C  c  2 ) 


(Ua) 


-  / 1  -  x 

/  v  - -  ^  C 

V  1  ♦  x, 

-  # 


1r  I  X  -*.l 


iNb) 


TT57 

- -  Jof  ( t 

1  *  t  J 


100 


<wi  f'J 


<  X 


rr'S't, 

~T~ 


A  '  2 


»n<3  .J. 


-X  )  - 


IX-X, 


d  X. 


dt,  5 


5  -  ex' 


Id 


and  tor  t .52  thev  are: 
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where  it  is  understood  that  ^  (i  - 1  )  <£0  *  [a’  -  O  c<rs  "  ’  x  a  J  .  It  so 

happens  that  both  the  radial  and  axial  integrations  defined  in  Eqns.  (13)  and 
(14)  can  be  performed  analytically,  which  is  a  great  advantage  over  schemes 
requiring  numerical  quadrature.  The  resulting  expressions  are  given  in  Appen¬ 
dices  B  and  C;  for  the  present,  it  is  assumed  that  these  quantities  are  all 
known. 


Next  we  choose  a  set  of  collocation  points  [  <Tn  ,  Xj  ]  ,  n  »  / ,  2.  ...NM 
Xs  1,2.. .ML  ,  at  which  to  evaluate  Eqn.  (11),  and  define 
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Substituting  Eqns.  (12)  -  (15)  into  Eqn.  (11)  reduces  the  integral  equation  to 

the  following  set  of  simultaneous  linear  algebraic  equations  for  the  CL  ■’  s  • 
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where  NN  *  NL  ■  N I  '  \J  t'or  the  system  to  he  determinant.  The  elements  of  the 
matrix  U  represent  the  dimensionless  upwash  produced  at  point  h\ 

the  (  l  y  a  )  loading  function.  These  are  given  bv : 


where  NM  and  NK  are  the  terms  at  which  the  m  and  k  series  are  truncated 


respectively 


In  writing  Eqn.  il"),  use  has  been  made  of  the  following  simplification 
in  the  wake  term.  For  each  i  i,  i  )»  it  becomes 


This  term  is  seen  to  vanish  identically  for  t he  uniform  spanwise  loading  function 
»  j  ,  as  it  should,  since  there  is  no  trailing  vorticitv  in  this  case. 


Before  solving  the  matrix  equation  in  (16),  one  more  feature  of  the 
solution  needs  to  be  considered;  namely,  the  hard-wall  boundary  condition  in 
Eqn.  (3)  requires  that  A-j have  :ero  spanwise  derivative  at  the  hub  and  tip. 

This  condition  can  be  enforced  automatically  as  follows.  If  Eqn.  (.12)  is 
differentiated  with  respect  to  cr»  »  and  the  resulting  expression  set  equal  to 
zero  at  <Ta*  -L  and  1,  two  equations  involving  the  will  follow.  Simultaneous 

satisfaction  of  these  two  equations  will ,  in  general,  depend  on  x „  .  But  if  it 
is  required  that  they  be  satisfied  for  all  l Xti  <1,  and  for  each  i  individually, 
then  (it  vj->  and  can  be  solved  for  uniquely: 

(13a) 

(13b) 

In  effect,  for  each  value  of  i  Eqn.  (5)  has  been  used  to  express  two  of  the 
radial  coefficients  in  terms  of  the  remaining  (  A /J-2). 
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Upon  substituting  Eqn.  (18)  into  (16),  one  obtains  an  altered  matrix 
equation  of  the  form: 


■VI  NS -1 


where  the  elements  of  the  modified  matrix  u' 


n  *  \,Z  .  .  .  NN 


are  given  by: 
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ML 


(19a) 


u- -  a;- f< r„  ,  x4) 


a  -•  r 


(ISb) 


U 


2  *  }  *  VJ  -i 


Note  that  the  ^  «i  terms  remain  unaffected,  since  they  satisfy  the  required  con¬ 
dition  trivially.  Since  the  number  of  unknowns  in  Eqn.  (19a)  has  effectively 
been  reduced  to  A/I  x  (IMJ  -  2  )  ,  the  number  of  collocation  points  should  be 
reduced  accordingly.  In  view  of  how  Eqn.  (IS)  was  derived,  the  most  logical 


approach  would  be  to  reduce  the  number  of  spanwise  stations  in  order  to  get 
NN  *ML  -  Ml  x(NJ-2)  .  Once  having  set  up  the  matrix  equation,  its  solution 
was  carried  out  by  using  Gaussian  elimination  with  pivoting  for  the  cases  reported 


here.  This  yields  CL  for  all  l  and  -j  -  '  •  •  •  M J  -2 

‘}  « 
are  then  obtained  from  Eqn.  (18). 


the  j-i  and  fl-i,  a/j 


The  local  loading  at  any  point  on  the  blade  is  calculated  from  Eqn. 

(12).  The  form  of  that  equation  also  allows  simple  analytical  expressions 
to  be  derived  for  the  sectional  lift  and  moment  coefficients,  the  overall  torque 
and  axial  force  coefficients,  the  static  and  total  pressure  rise,  and  the  turning 
produced  by  the  rotor.  The  integrations  involved  are  elementary,  and  so  only 
the  results  will  be  quoted. 


The  sectional  lift  on  the  blade  Lear')  is  perpendicular  to  the  local 
undisturbed  flow,  and  considered  positive  when  in  the  tn  direction  (Fig.  2). 
We  define  the  lift  coefficient  as 


(<r)  s 


L  c<r; 

4-  /»—  u  £*. 


0  +(4>r  cr)22Vz 


I  a-o  C<r,  .*)  a  x 


TT  [i +(0T  <r)*]’/z  J. 

- 2 -  H, <r 


T*.,> 


The  pitching  moment  M(<r)  is  taken  about  the  blade  mid-chord  and  considered 
positive  when  clockwise  (Fig.  2).  Its  coefficient  is  defined  as 


■Just  as  in  isolated  airfoil  theory,  wo  see  that  only  the  first  few  chordwise 
load  functions  contribute  to  the  net  lift  and  moment. 


The  overall  torque  exerted  on  the  fluid  by  the  complete  blade  row,  T 
is  considered  positive  when  counter-clockwise  as  one  looks  downstream  (Fig.  1). 
he  define  the  torque  coefficient  as: 


Z  4.  r 


f  irCLcr)cL<r 

J,  [’  +  Cv>,  r)  *]  ''2 


-  f  i1  >(is±l) 

l •»  *  J  l  j  *  /  / 

The  axial  farce  exerted  by  the  blade  ro»  on  the  fluid,  F*  U  considered 
positive  "hen  it  points  upstream.  The  analogous  coefficient  is  obtained  from: 


C  T  CLi<r)  a  t 
[i  <r)  tJ  ’'*■ 


The  net  static  pressure  rise  across  the  row  can  he  evalauted  from 

Eqn*  U)  at  1  *  *  ®  •  and  when  normalized  by  the  dynamic  pressure  of  the  axial 

flow  becomes 
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The  static  pressure  ratio  is  easilv  shown  to  be- 

where  here  3 ’  is  the  specific  heat  ratio. 

In  duct  coordinates,  the  rise  in  total  enthalpy,  and  hence  total 
temperature,  can  be  related  to  the  work  done  on  the  fluid.  The  total  pressure 
ratio  can  then  be  obtained  from  the  isentropic  relations,  and  the  result  is. 
to  first  order  in  small  quantities: 


<*«„>  M  1  _  *  3  Pr  n-  ^3CL(a-)  c*3 


1°*.c 
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The  brackets  on  the  left  indicate  that  a  circumferential  average  has  been 
taken.  After  allowing  for  differences  in  notation,  it  is  readily  shown  that 
Eqns.  (24)  and  (26)  are  equivalent  to  those  given  in  Ref.  3. 

One  can  also  express  the  circumferentially-averaged  turning  produced 
by  the  rotor  in  either  duct  or  blade-fixed  coordinates,  or  respectively, 
in  terms  of  the  above  parameters.  Turning  is  considered  positive  when  the  flow 
is  deflected  counter-clockwise  in  Fig.  2.  To  first  order,  AD--<.'irs>/U  and 
AD  =  -  <'J„> /  Uh  ,  where  V&  can  be  put  in  terms  of  and  -p  through: 

=  irn  yr  +  zrs  yr 

-  0  +  <P*  <r*)  'z  [  vn  +  (<pr  <r)  v5 ]  ( 27) 

■P 

^  i-8> 

and  -p  and  trn  are  given  by  Eqns.  (1)  and  (4),  respectively.  The  results  are: 


8  pr  CLC<r) _  Hip 

2  ir[j  +(<fir  r)*]’'*  “  UL9t<T) 

nr  1*  <r  cy _ 

0 


(291 

(30) 


It  is  seen  that,  whereas  A0  depends  only  on  conditions  at  the  local  radius, 
-ig  is  influenced  by  the  overall  torque  as  well.  Eqn.  (50)  is  equivalent  to 

3 

that  cited  by  McCune  and  Dharwadkar. 


D.  NUMERICAL  RESULTS  AND  DISCUSSION 

A  computer  program  implementing  the  solution  procedure  described  above 
has  been  developed  using  double  precision  arithmetic  throughout  so  as  to 
minimite  the  influence  of  round-off  error.  Six  cases  have  been  run  to  test 
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the  program  over  a  broad  range  of  operating  conditions,  as  specified  by  the 
input  parameters  B  ,  ,  M  ,  (pT  and  iq  r  .  Table  1  summarizes  the  values 

used  for  each  case,  along  with  a  brief  description  of  the  specified  camber 
line.  The  first  three  cases  were  designed  to  test  the  present  theory  against 
a  two-dimensional  strip  theory  at  a  high  hub/tip  ratio,  where  one  would  expect 
strip  theory  to  provide  a  reasonable  approximation.  A  similar  comparison  is 
drawn  in  Case  4  for  a  low  hub/tip  ratio,  for  which  three-dimensional  effects 
are  much  more  significant.  Finally,  Cases  5’ and  6  provide  comparisons  with 
the  inverse  three-dimensional  theory  of  Okurounmu  and  McCune.^’7 

TABLE  1 

INPUT  PARAMETERS  FOR  CASES  1-6 


Case  * 

3 

i.  Si.  1 

'  ~  i 

Axial  Mach  * 

.  ujrr 

O  x  —  ■ 

■r  (J 

r 

It  2  rr 

j 

I 

Camber  Line  Shape  | 

! 

1 

Ut 

Cj 

0.3 

Q.S 

1.0 

0.03 

Flat  plate  at  5° 
local  angle  of  attack ■ 

2 

50 

0.3 

0.3 

1.0 

0.03 

10*.  parabolic-arc 
camber 

3 

j 

50 

0.3 

0.5 

1.0 

0.06 

Fiat  plate  at  3a 
local  angle  of  attack j 

i  4 

1 

30 

0.3 

0.5 

1.0 

o.os 

Flat  plate  at  3* 
local  angle  of  attack1 

3 

1 

40 

0.3 

C  .5 

1.497 

0.0833 

Free  vortex  design 

j  6 

40 

0.3 

O.S 

1 .497 

0  .0833 

20%  variation  in  P 

In  each  case,  NI  =  5  chordwise  and  N'J  =  5  spanwise  loading  functions 
were  used  to  represent  A-p  .  Along  the  chord,  NL  =  S  collocation  stations  were 
equally-spaced  between  the  leading  and  trailing  edges,  X  =  il;  in  the  spanwise 
direction,  NN  =  3  collocation  stations  were  equally-spaced  between  the  hub  and 
tip,  <r  =  h,  1.  The  number  of  azimuthal  modes  included  in  the  calculation, 

MM,  was  10  in  all  cases;  the  number  of  radial  modes  included,  NK,  was  set  at 
20  or  30  depending  on  whether  A  was  0.S  or  0.5,  respectively. 


1. 


Comparisons  With  2-D  Strip  Theory 


In  Case  I,  the  camber  line  was  specified  to  be  a  flat  plate  at 

each  radial  station;  however,  a  radial  twist  is  imparted  to  the  blades  such  as 

to  hold  the  local  angle  of  attack,  a J0  ,  constant  at  *5  degrees.  The  high 

hub/tip  ratio  in  this  case  suggests  that  the  loading  might  be  reasonably  well 

predicted  by  application  of  a  two-dimensional  strip  theory.  However,  caution 

must  be  exercised  in  making  such  comparisons;  as  pointed  out  by  McCune  and 
3 

Dharwadkar,  two-dimensional  cascade  theory  defines  the  "undisturbed”  flow- 
velocity  as  the  vector  mean  of  the  inflow  and  outflow  velocities,  whereas  in 
the  present  three-dimensional  theory,  it  is  defined  to  be  the  inflow  velocity, 

LL  .  Since  the  difference  between  these  reference  directions  does  not 
depend  on  axial  position,  it  amounts  to  using  a  different  angle  of  attack  in 
the  two-dimensional  theory,  a,0  .  than  would  be  used  in  the  three-dimensional 
theory,  x30  ,  though  the  physical  blade  geometry  and  inflow  velocity  are  the 
same.  The  difference  between  the  two  angles  is  proportional  to  the  net  turning 
across  the  row,  which  in  turn  can  be  related  to  the  sectional  loading: 

L  (.31) 

r  ,  .  J~ 

— r~m  w 

The  lift  coefficient  CL  is  more  convenient  to  work  with  in  a  strip  theory  than 
that  defined  in  Eqn.  (20).  It  is  easily  shown  that  to  first  order  in  the  per¬ 
turbation  scheme,  the  two  are  related  by  CL-  CL  cos'  V  . 


’20 


’30 


3; 

4  rr  cr  3  * 


Now  the  two-dimensional  theory  predicts  that 


‘  -3,  U  ’ 

where  the  factor  is  defined  as  the  ratio  of  the  lift  on  a  blade  in  cascade 
to  the  lift  on  the  same  blade  operated  as  an  isolated  airfoil.  For  flat  plate 

blades  in  incompressible  flow,  this  factor  has  been  expressed  in  closed  form  b\ 

•  17*  c 

Pistolesi  as  a  function  of  the  stagger  angle,  o  ,  and  solidity,  ; 
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The  values  of  stagger  and  solidity  to  be  used  in  Eqn.  i35)  are  those  in  the 

corresponding  incompressible  flow.  They  are  related  to  the  values  in  the 

original  cascade  flow  through  the  Prandt 1 -Glauert  transformation  described  b\ 
13 
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Elimination  of  xi0  between  Eqns.  (31)  and  (52)  then  gives: 
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A  comparison  of  the  strip  theory  predictions  of  F.qn.  ioSI  for  the 
conditions  of  Case  1  with  our  three-dimensional  calculations  is  given  in  Fig.  5. 
When  viewing  these  and  subsequent  results  for  CL  it  should  be  kept  in  mind  that 
the  spanwise  variation  in  L  will  be  even  greater,  since  the  local  dynamic 
pressure  and  chord  both  increase  with  radius.  The  agreement  between  the  two 
calculations  is  seen  to  be  quite  good  at  mid-annulus.  However,  whereas  the 
strip  theory  predicts  a  5 Si  increase  in  •CL  from  hub  to  tip,  the  three-dimensional 
results  show  little  variation.  The  corresponding  changes  in  the  circulation,/"'  , 
from  hub  to  tip  are  +6"%  and  19%,  respectively.  Thus,  significant  trailing 
vortex  effects  can  be  expected  despite  the  high  hub 'tip  ratio. 

The  cascade  factor  in  Eqn.  (35)  was  computed  to  be  0.9b  at  mid¬ 
annulus,  indicating  tthat  the  blades  are  operating  nearly  as  isolated  airfoils. 
This  results  from  the  low  solidity  of  the  blade  row,  which  at  this  radius  is 
roughly  two-thirds.  Hence,  one  would  expect  the  chordwise  load  distribution 
to  closely  approximate  that  of  an  isolated  flat  plate.  This  comparison  is  shown 
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L 


0. 


in  Fig.  4.  where  the  isolated  airfoil  result  has  been  normal i:ed  b>  the  factor 
so  as  to  yield  the  same  sectional  lift,  i.e.. 


y  p-*>  U* 


•4  01  >  p  .  P 

-  uv**  — 

3,  ^ 


The  agreement  is  seen  to  be  excellent 


In  Case  2,  the  blades  are  assumed  to  have  :ero  angle  of  attack  relative 
to  the  inflow  velocity  and  a  parabolic-arc  camber  line;  otherwise,  r ho  input 
parameters  are  the  same  as  in  Case  1.  The  ratio  of  maximum  camber,  ,  to 

the  local  chord  length  is  take  to  be  0.10,  and  constant  with  radius.  \s  m 
Case  1,  one  would  expect  that  the  high  hub  tip  ratio  and  low  solidity  would 
cause  the  three-dimensional  results  to  approach  those  for  a  two-dimensional 
isolated  airfoil.  In  the  two-dimensional  theory,  the  blade  camber  contributes 
to  the  lift  coefficient  in  the  amount 

-  2**0. 

L,  * - - 

p* 

where  is  the  cascade  correction  factor  for  camber  analogous  to  in  the 

case  of  the  flat  plate  cascade.  Unfortunately,  while  kA  can  be  expressed  in 
closed  form,  K  ^  cannot.  But  it  was  seen  above  that  k  ^  is  near  unit)  anvwn>  , 
and  so  it  seems  reasonable  to  assume  that  for  these  conditions  )iu  is  also.  It 
should  also  be  noted  that  even  though  ^  =  0  in  this  case,  remains  finite 

since  there  is  still  a  finite  turning.  Thus,  if  the  contributions  from  tqns.  (,521 
and  (.3?)  are  added,  and  Eqn.  (.'ll  used  to  eliminate  ol£D.  the  expression  that 
results  for  Cu  is 


The  comparison  between  the  lift  coefficients  predicted  by  Eon.  (58)  and 
the  full  three-dimensional  theory  for  Case  2  is  also  shown  in  Fig.  5.  Vgain, 
good  agreement  is  shown  at  mid-annulus,  while  the  strip  theory  significantly 
overestimates  the  lift  near  the  tip,  and  underestimates  it  near  the  hub.  Having 
obtained  the  two-dimensional  Cu  from  Eqn.  (53),  <xeo  can  be  determined  from 
Eqn.  (51).  One  can  then  compute  the  corresponding  chordwise  load  distribution, 
which  is  a  superposition  of  the  contributions  from  angle  of  attack  and  camber. 
The  comparison  between  this  prediction, 
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and  the  three-dimensional  results  at  mid-annulus  is  shown  in  Fig.  4,  and  again 
very  good  agreement  is  obtained. 


For  Case  5  ,  and  hence  the  solidity  all  across  the  annulus,  has 

been  doubled;  the  other  input  parameters  and  the  camber  line  are  the  same  as  in 
Case  1.  This  case  thus  provides  a  better  check  on  whether  the  present  program 
correctly  predicts  blade  interference  effects.  Its  predictions  and  those  from 
the  strip  theory,  Eqn.  (55),  are  also  compared  in  Fig.  5.  The  higher  solidity 
is  seen  to  result  in  an  overall  reduction  in  loading,  the  factor  now  being 

only  0.65  at  mid-annulus.  Again  very  good  agreement  between  the  two  theories 
is  exhibited  at  mid-annulus,  with  the  discrepancies  at  other  stations  following 
the  same  trends  seen  in  Cases  1  and  2.  Whereas  in  those  cases,  the  chordwise 
load  distributions  closely  approximated  those  for  an  isolated  airfoil,  in  this 
case  one  would  expect  the  higher  solidity  (about  1.5)  to  alter  the  distribution 
somewhat.  This  is  illustrated  in  Fig.  5  where  the  present  results  are  again 
compared  with  Eqn.  (56).  Figs.  5  and  5  together  show  that  the  increased  solidit 
acts  to  reduce  the  overall  lift,  and  to  shift  more  of  the  loading  forward  near 
the  leading  edge  at  the  expense  of  that  over  the  rear  half  of  the  blade. 


The  throe  cases  discussed  thus  far  were  run  to  assess  the  agreement 
between  the  present  analysis  and  strip  theory  for  high  hub/tip  ratio.  The  con¬ 
ditions  for  Case  4  were  chosen  to  exhibit  more  significant  three-dimensional 
effects,  with  the  hub  tip  ratio  reduced  to  0.5  and  the  number  of  blades  reduced 
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to  30;  the  other  parameters  and  the  camber  line  are  the  same  as  in  Case  1.  The 
results  and  comparison  with  strip  theory  are  shown  in  Fig.  b.  Both  calculation 
show  larger  radial  variations  than  in  the  previous  cases,  as  expected,  \gain, 
the  strip  theory  significantly  underestimates  the  loading  at  inboard  stations, 
and  overestimates  it  near  the  tip,  with  good  agreement  at  mid-annulus.  However 
the  percentage  errors  relative  to  the  three-dimensional  predictions  are  much 
greater  in  this  case  than  in  the  previous  cases  with  high  hub/tip  ratio.  It 
should  be  noted  here  that  the  strip  theory  predictions  from  cr  -  0.8  outward 
are  identical  to  those  for  Case  1.  This  occurs  because  the  two  parameters  3 
and  nr  enter  the  strip  theory  calculation  only  through  the  cascade  solidity, 
and  then  only  as  the  product  3  rj r  ,  which  is  the  same  in  both  cases. 

The  local  chordwise  load  distributions  at  the  hub,  mid-annulus,  and 
tip  radii  as  predicted  by  the  three-dimensional  theory  are  compared  in  Fig. 
Surprisingl v,  there  is  seen  to  be  very  little  variation  between  the  three.  \'o 
doubt  this  is  due  in  large  part  to  our  having  twisted  the  blades  in  such  a  way 
as  to  maintain  the  geometric  angle  of  attack  constant  with  radius. 

In  all  four  cases  discussed  thus  far,  it  is  seen  that  the  present 

three-dimensional  theory  and  a  two-dimensional  strip  theory  show  good  agreement 

only  near  mid-annulus.  The  strip  theory  tends  to  underestimate  the  loading  neai 

the  hub  and  overestimate  it  near  the  tip,  even  when  the  hub/tip  ratio  is  fairly 

§ 

high.  McCune  and  Dharwadkar'  have  found  the  same  effect  in  their  study  of  the 
three-dimensional  lifting-line  approximation.  In  their  model,  each  blade  was 
represented  by  a  single  radially-oriented  bound  vortex.  The  strength  of  the 
vortex  was  allowed  to  vary  radially,  thus  generating  a  trailing  vortex  pattern. 
By  computing  the  resultant  induced  angle  of  attack  at  each  spanwise  station, 
McCune  and  Pharwadkar  were  thus  able  to  correct  the  sectional  lift  predicted 
by  strip  theory  to  account  for  this  three-dimensional  effect.  They  concluded 
that  the  distribution  of  trailing  vorticitv  acts  in  such  a  way  as  to  reduce 
spanwise  variations,  which  is  consistent  with  the  results  presented  here  in 
Figs.  5  and  6. 
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2 .  Comparisons  With  Inverse  3-D  Theory 

The  remaining  two  cases,  5  and  o,  are  of  particular  interest 
because  they  provide  comparisons  between  the  present  work  and  another  three- 
dimensional  lifting-surface  theory,  that  of  Okurounmu  and  McCune.0’  In  contrast 
to  the  off-design  problem  considered  here,  they  solved  the  design  problem,  i.e., 
the  determination  of  that  camber  line  which  will  produce  a  specified  load  dis¬ 
tribution.  They  presented  plots  of  camber  line  slope  vs.  chordwise  distance 
which  we  input  to  our  program  to  see  if  the  predicted  load  agreed  with  that 
specified  in  the  design  calculation.  Two  calculations  were  performed;  the  first 
was  for  a  free-vortex  design  (Case  5),  in  which  by  definition  the  circulation 
about  the  blades,  P*  ,  was  constant  with  radius  (Okurounmu  and  McCune  refer 
to  this  as  the  constant-work  case'.  In  the  second  (Case  6),  the  circulation 
varied  by  20*  from  hub  to  tip,  referenced  to  its  radius-weighted  mean  over  the 
span.  In  both  cases,  the  variation  in  design  total  pressure  ratio  across  the 
annulus  was  in  the  range  1.10  to  1.15.  The  camber  line  slopes  had  been  normalized 
by  the  overall  static  pressure  rise  coefficient  across  the  rotors,  as  defined 
by  Eqn.  (24).  This  quantity  was  assumed  to  be  unity  in  reading  the  curies,  so 
that  the  present  program  should  predict  a  value  of  of  one  if  the  two 

theories  are  truly  consistent. 

For  Case  5,  the  present  program  predicts  a  value  of  =  1.03', 

which  is  considered  very  good  in  view  of  the  accuracy  with  which  the  camber 
line  slopes  could  be  read  from  the  graphs.  T’he  predicted  circulation,  norma¬ 
lized  by  the  product  of  the  constant  axial  velocity,  U  ,  and  the  tangential 
blade  spacing  at  the  tip,  Lq  ,  is  compared  with  the  constant  design  value  in 
Fig.  3.  Good  agreement  is  exhibited  over  the  entire  span.  The  distributions 
of  CL  are  also  compared,  the  design  values  having  been  inferred  from  Eqn.  (29). 
This  quantity,  rather  than  CL  ,  was  plotted  because  it  reflects  the  true  radial 
variation  in  lift. 

Fig.  9  compares  the  chordwise  load  distribution  at  mid-annulus  as 
predicted  by  the  present  analysis  with  the  specified  design  distribution.  The 
latter  was  proportional  to  (i  -x2>  e  at  each  radial  station.^'  Good  agreement 
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is  shown  over  most  of  the  chord;  in  principle,  one  should  be  able  to  reduce 
the  discrepancies  near  the  leading  edge  by  increasing  the  number  of  chordwise 
collocation  points.  However,  such  a  calculation  was  not  felt  warranted  in 
view  of  the  inherent  inaccuracy  in  reading  the  slopes,  as  noted  above. 

Unfortunately ,  the  results  for  Case  6  do  not  exhibit  good  igreement. 
In  this  case,  was  predicted  to  be  1.34,  or  more  than  50^  too  high.  This  is 

reflected  in  Fig.  10,  which  compares  the  predicted  spanwise  variations  in 
and  CL  with  those  prescribed  by  Okurounmu  and  McCune;  our  results  are  signifi¬ 
cantly  higher  than  the  design  distributions  over  the  entire  span. 

The  disagreement  in  this  case  has  not  yet  been  explained;  it  is  too 
large  to  be  reasonably  ascribed  to  errors  in  reading  the  camber  line  slopes.  As 
seen  from  Table  1,  the  inputs  for  Cases  5  and  b  are  identical  except  for  the 
blade  geometry.  In  Case  6,  this  leads  to  significant  trailing  vortex  effects 
not  present  in  Case  5.  The  vortex  representation  of  Refs.  6  and  ~  uses  an 
entirely  different  analytical  expresssion  for  the  wake  terms  than  the  dipole 
representation  used  here,  as  discussed  in  Appendix  A.  However,  trailing  vortex 
effects  were  by  definition  excluded  a  priori  from  the  free-vortex  design  of 
Refs,  b  and  7.  Under  these  circumstances ,  the  most  logical  explanation  for 
the  good  agreement  in  Case  5  and  lack  of  agreement  in  Case  o  would  be  an  error 
in  the  numerical  evaluation  of  the  wake  terms. 

In  defense  of  the  present  calculations ,  it  should  be  noted  that  the 
same  number  of  se  loading  functions  was  used  to  represent  Ap  in  both 

Cases  5  and  6;  in  fact,  the  matrix  equations  which  were  solved  in  each  case, 

Eqn.  (19),  were  identical  except  for  the  column  vector  of  known  slopes  on  the 
right-hand  side.  Thus,  trailing  vorticity  had  the  same  opportunity  to  influence 
the  results  of  Case  5  as  Case  6,  although  in  the  former  it  correctly  proved 
unimportant.  If  an  error  were  present  in  our  evaluation  of  these  terms,  it  is 
difficult  to  see  how  it  could  produce  significant  discrepancies  in  one  case  and 
not  the  other.  It  is  hoped  to  shed  further  light  on  this  question  in  the  future. 
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An  investigation  has  been  carried  out  on  the  effect  of  vary  ini; 
the  number  of  itimuthal  and  radial  modes  included  in  the  calculation ,  ,\jm  and 
NK  .  The  static  pressure  rise  coefficient  defined  in  Iain,  (.24J  was  chosen  as 
a  representative  indicator  of  the  overall  convergence  of  the  calculations.  The 
convergence  of  with  A/M  is  shown  in  t  ig.  11  for  Cases  l  and  I.  The  value 

at  MM  =  10  was  chosen  as  a  reference,  and  NK  was  held  constant  at  a  value  of 
20  for  Case  1  and  50  for  Case  I .  It  is  seen  that  the  convergence  is  fairly 
rapid,  and  that  as  few  .is  4  atimuthal  modes  give  reasonably  accurate  results. 
Also,  the  greater  three-dimensional itv  in  Case  J  appears  to  aid  the  convergence 
in  Fig.  11,  but  whether  this  would  hold  true  in  general  is  hard  to  say.  Fig.  12 
shows  the  convergence  for  Case  5  with  NK  =  20;  here  the  actual  values  of  Cps 
are  plotted  because  the  relative  variations  arc  much  greater,  and  the  "correct" 
or  design  value  was  known  a  priori.  It  appears  that  in  this  case,  which  has  a 
solidity  more  than  double  that  of  Cases  1  or  4,  at  least  9  or  10  terms  are 
needed  to  get  the  error  in  below  1 . 

The  convergence  of  the  calculations  with  NK  is  shown  in  Fig.  15 
for  Cases  1,  4  and  5.  The  value  of  MM  was  held  constant  at  1 0  in  all  runs. 
Noting  the  different  scale  factors  on  each  curve,  it  is  seen  that  for  comparable 
accuracy  one  must  include  many  more  radial  modes  in  the  highly  three-dimensional 
Case  4  than  the  more  two-dimensional  Case  l,  which  is  intuitively  satisfying, 
lhe  highest  solidity  case,  Case  5,  again  exhibits  the  slowest  convergence. 

It  should  be  kept  in  mi  ml  that  the  convergence  of  more  localized 
quantities,  e.g.,  the  sectional  lift,  L  ,  or  the  local  load  distribution,  .A  p  , 
can  lie  expected  to  be  slower  than  that  exhibited  for  .  For  Cases  1,  4  and 

S,  these  appear  to  have  converged  within  1A  with  respect  to  NK  .  Cases  1  and  I 
hail  also  converged  to  this  extent  with  respect  to  NM  ,  although  Case  5  still 
exhibited  variations  on  t he  order  of  S-9A  between  MM  =  9  and  10.  At  present, 
the  program  storage  is  limited  to  NM  $  10,  and  its  modification  to  allow  further 
increases  in  MM  for  this  one  case  does  not  seem  warranted  at  this  time. 


In  summary,  it  appears  that  the  inclusion  ot‘  10  azimuthal  modes  cut 
between  20  ind  30  radial  modes,  depending  on  hub  t  i  ]'  ratio,  provides  sufficient 
accuracy  tor  most  purposes,  at  least  with  regard  to  measures  of  overall  perform 
ance.  Of  course,  this  conclusion  may  have  to  be  modified  for  parameters  signifi 
cantly  outside  the  range  covered  in  Table  1.  Qualitatively,  it  appears  that 
high  solidity  aggravates  the  convergence  of  both  the  azimuthal  and  radial  mode 
sen's.  low  hub  tip  ratio  cases  require  more  radial  modes,  as  expected,  but 
may  hate  a  small  advantageous  effect  on  the  convergence  of  the  azimuthal  mode 
series. 

The  computer  running  times  required  for  Cases  l-o  varied  between 
10  and  20  minutes  on  m  IBM  5o0,  o5,  using  the  G  compiler.  It  is  estimated  that 
the  running  time  on  an  IBM  370  would  be  reduced  b>  a  factor  of  6  or  ”,  U>  far 
the  maiorit>  of  time  is  spent  in  evaluating  the  axial  integrals.  The  time 
needed  for  their  evaluation  is  primarily  dictated  by  the  number  of  integrals 
which  can  be  done  analytically  before  round-off  error  becomes  a  problem  (see 
Appendix  Cl.  \s  a  result,  the  running  time  can  be  expected  to  increase  along 
with  the  solidity.  Vs  noted  in  connection  with  Cases  5  and  o,  however,  the 
coefficient  matrix  in  tqn.  (I'A1  depends  only  on  the  first  five  parameters  in 
Table  l.  The  program  lias  provision  for  storing  the  matrix  on  tape,  so  that 
subsequent  calculations  for  the  same  operating  conditions,  but  different  blade 
camber  lines,  can  be  performed  in  a  matter  of  seconds. 


i: .  CONCLUSIONS 

V  method  has  been  presented  for  computing  the  stead'  loading  on  a 
three-dimensional  annular  Made  row  in  compressible  flow,  given  the  blade 
camber  line.  V  kerne l -function  procedure  is  used  to  solve  the  integral  equation 
relating  the  unknown  loading  to  the  camber  line  slope.  The  greatest  advantage 
of  the  solution  procedure  presented  here  is  that  if  allows  analvtical  expressions 
to  he  derived  tor  both  tlve  spanwise  and  chordwiso  integrals  which  are  required. 

It  also  allows  the  various  aerodynamic  and  performance  parameters  of  t  lie  rotor 
to  be  expressed  as  simple  algebraic  functions  of  the  loading  expansion  ce 
e ft'  i  c  i  ent  s  . 


Calculations  made  thus  tar  have  concentrated  on  comparing  the  present 
predictions  against  those  ot'  two-dimensional  s  t  rip  theory,  is  well  as  the  three 
dimensional  inverse  1  iftmg-surt.ice  theory  of  Okurounmu  and  Mciamo.  '  dood 
agreement  between  t ho  three-dimensional  and  strip  theories  occurs  only  near 
mid-annulus,  even  for  a  moderately  high  hub  tip  ratio.  For  low-solidit>  blade 
rows,  it  is  also  demonstrated  that  the  present  analysis  predicts  chordwise 
loading  distributions  at  mid-annulus  which  are  in  excellent  agreement  with  those 
of  isolated  thin-airfoil  theory,  for  blades  with  both  angle  of  attack  and  camber. 
The  three-dimensional  theory  predicts  loadings  which  ire  significant lv  higher 
near  the  hub  and  lower  near  the  tip  than  the  strip  theory,  thus  giving  much 
lower  radial  variations.  Such  behavior  can  be  attributed  to  the  induced  angle 
of  attack  distribution  generated  by  the  trailing  vortex  wakes,  which  is  not 
accounted  for  in  strip  theory;  this  smoothing  appears  to  be  the  principle  effect 
of  three-dimensional ity,  at  least  for  subsonic  relative  Mach  numbers. 

Comparison  with  the  inverse  three-dimensional  theory  of  Okurounmu  and 
McCune*1’  shows  good  agreement  over  the  whole  span  for  the  case  in  which  the 
circulation  is  constant  with  radius.  The  chordwise  load  distribution  at  mid¬ 
annulus  also  compared  well  with  the  variation  prescribed  by  them.  However,  for 
blades  designed  to  produce  a  linear  variation  in  circulation  from  hub  to  tip, 
the  present  theor>  predicts  loadings  which  are  s ign i f icant l>  higher  than  the 
design  distribution.  It  is  hoped  to  resolve  this  discrepancy  in  the  future. 

When  coupled  with  the  analogous  treatment  of  the  disturbance  flowfield 
2  -  4  1J 

produced  by  blade  thickness,  '  the  loading  analysis  presented  lie  re  should 
provide  an  efficient  means  of  studying  three-dimensional  effects  in  turbo 
machinery  rows.  I'he  program  should  prove  particularly  advantageous  for  evaluating 
alternate  blade  camber  lines  at  a  given  operating  condition.  The  coefficient 
matrix  in  Fqn.  (l'.t-)  then  needs  to  be  evaluated  and  inverted  onl>  once,  so  that 
results  for  additional  blade  camber  lines  can  bo  obtained  in  a  matter  of  seconds, 
whereas  a  finite-difference  solution,  for  example,  would  require  a  whole  new 
calculation  for  each  case.  Thus,  the  program  could  serve  to  narrow  down  the 
number  of  candidate  blade  profiles  before  going  to  a  f i nite-d i f t'erence  solution 
of  the  full  nonlinear  equations.  It  could  also  prove  useful  in  providing 


better  initial  conditions  to  a  difference  scheme.  illowing  t ho  latter  to  converge 
in  fewer  iterations.  Some  development  work  remains,  however.  In  particular, 
the  effect  of  varying  five  position  nd  number  of  the  collocation  :vmt<,  and 
hence  the  number  of  terms  in  the  expansion  of  t lie  unknown  loading,  should  be 
investigated.  \lso,  the  assumption  of  constant  axial  chord  projection  should 
be  relaxed  to  allow  application  to  a  wider  varietv  of  blade  geometries. 

Initially,  it  was  planned  to  extend  the  reported  1  i  ft  mg- surf 'ace  m.ilvsis 
to  supersonic  tip  speeds.  However,  it  has  been  decided  instead  to  proceed  with 
the  extension  to  unsteadv  flow  it  subsonic  tip  speeds,  for  the  following 
reasons.  In  previous  linear; ted  treatments  of  both  the  thickness"  1  and  inverse- 

loading''  problems,  a  transonic  resonance  has  been  shown  to  occur  it  super¬ 

sonic  tip  speeds.  The  amplitude  of  each  duct  acoustic  mode  is  predicted  to 
become  infinite  at  its  cut  off  condition,  though  in  reality  the  effects  of 
viscosity  and  nonlinearities  in  the  flow  would  prevent  this  from  happening. 

\ttempts  have  been  made  to  incorporate  viscous--1  ,md  nonlinear  effects’^  into  the 
linearized  model  in  an  approximate,  ad  hoc  fashion.  Uthough  the  resonant  modes 
can  be  successfully  limited  in  this  wav,  their  resultant  amplitudes  are  ten  much 
dependent  on  the  particular  assumptions  made  concerning  the  damping  and  or  non¬ 
linear  effects. 

\  more  rigorous  approach  to  the  problem  at  supersonic  tip  speeds 
would  retain  the  nonlinearities  in  the  governing  equations  to  be  solved,  rather 
than  trying  to  approximate  them  in  a  linearized  analysis.  This  of  course 
requires  numerical  solution  procedures,  such  as  finite  difference  techniques, 
rhree-dimcnsional  relaxation  solutions  of  the  transonic  nonlinear  small  di'v 

turbance  equations,  which  are  capable  of  including  weak  shocks  in  t he*  flow,  have 

.  I  o 

been  obtained  tor  subsonic,  but  supercritical,  flow  through  a  rotor.  tfforts 

are  currently  underway  to  extend  this  technique  to  numerical  computations  for 

supersonic  relative  intlows.  Such  an  approach  not  only  avoids  the  problem  of 

transonic  resonance,  but  is  also  extendable  to  the  full  nonlinear  equations 

of  motion,  including  the  possibility  of  strong  shock  waves.  \  computer  code 

for  solving  the  full  three-dimensional  Tailor  equations  for  the  stead'  flow 

through  a  rotor  is  presently  being  developed  at  Cal  span  under  another  program."' 


For  the  above  reasons,  it  was  felt  that  a  more  useful  extension  of  the 
linearized  three-dimensional  analysis  would  be  to  the  study  of  subsonic  unsteady 
flow  through  a  rotor.  The  use  of  small -perturbation  analyses  in  the  study  of 
the  blades'  aeroacoustic  response  to  inflow  distortion  and  to  the  problem  of 
blade  flutter  is  well  established.  This  effort  is  described  in  Section  II. 
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SECTION  II 


UNSTEADY  LIFTING -SURFACE  THEORY  FOR  AN  .ANNULAR  ROTOR 

A.  INTRODUCTION 

The  previous  section  dealt  with  the  prediction  of  the  steady  loading 
on  an  annular  rotor  for  which  the  blade  camber  line  has  been  specified.  In 
this  problem,  the  flow  is  viewed  as  stationary  in  blade-fixed  coordinates  on 
the  assumption  that  the  rotor  is  isolated,  is  operating  in  a  uniform  inflow, 
and  consists  of  perfectly  rigid  blades.  In  practice  of  course,  such  conditions 
are  never  perfectly  realized.  The  close  proximity  of  other  blade  rows,  as  well 
as  non-uniformities  in  the  incoming  air,  give  rise  to  distorted  inflows  which 
generate  unsteady  blade  forces.  These  in  turn  can  lead  to  degradations  in  per¬ 
formance.  Also,  vibrations  of  the  blades  themselves  can  sometimes  couple  to 
the  resulting  unsteady  aerodynamic  field  in  such  a  way  as  to  lead  to  serious 
aeroelastic  problems.  Unsteady  aerodynamic  loads,  whether  from  inflow  dis¬ 
tortion  or  blade  flutter,  also  act  as  efficient  acoustic  sources,  and  hence 
their  reduction  is  one  means  of  reducing  engine  noise. 

With  the  above  applications  in  mind,  work  was  begun  on  the  formula¬ 
tion  of  an  unsteady  three-dimensional  lifting-surface  theory  for  an  annular 

1  O 

Made  row.  Namba'  also  has  extended  tyis  analysis  of  the  stead)  flow  case  to 

treat  the  unsteady  problem,  kobavashi  later  applied  Namba ' s  analysis  to  the 

study  of  pure  tone  fan  noise  due  to  inflow  distortion.  However,  the  same  scale 

factor  appears  in  this  formulation  as  was  included  in  the  steady  problem.  As 

mentioned  in  Section  l-B,  we  have  found  no  justification  for  its  presence,  and 

hence  feel  the  numerical  results  in  Refs.  JS  and  29  must  be  viewed  with  caution 
1  •>  13 

Salaun  " ’  has  also  presented  a  lifting-surface  theory  for  the  unsteady  case, 
as  noted  earlier,  and  we  have  shown  that  our  expression  for  the  pressure  field 
generated  by  the  blade  row  is  in  agreement  with  his.  We  are  hopeful  that  the 
application  of  the  techniques  developed  earlier  for  the  steady  problem  will 
lead  to  an  efficient  inversion  of  the  unsteady  integral  equation  as  well. 


Section  1 1  - B  below  presents  the  assumptions  made  in  the  analysis 
and  the  governing  partial  differential  equation.  In  Section  1I-C,  solutions 
to  the  equation  are  given  corresponding  to  unsteady  pressure  monopole  (.source! 
and  dipole  (force!  singularities.  The  latter  are  used  in  Section  II-P  to 
develop  the  pressure  field  for  the  entire  blade  row.  This  is  used  in 
Section  II-F.  to  derive  the  governing  integral  equation,  whose  features  are 
then  discussed  relative  to  the  proposed  solution  procedure.  Wherever  possible, 
the  notation  has  been  kept  consistent  with  that  used  in  Section  1  for  the 
steady  case,  with  a  few  exceptions  which  will  be  noted  as  they  arise. 


B.  FLOW  NR3PEL  AND  GOVERNING  EQUATIONS 

The  same  geometry*  is  assumed  for  the  present  case  as  had  been  used 
in  our  earlier  study  of  the  steady  load  problem.  That  is.  the  blade  row  is 
assumed  to  be  housed  in  an  infinitely  long,  hard-walled  annular  duct  of  con¬ 
stant  hub/tip  ratio,  h  ,  containing  a  uniform,  inviseid,  subsonic  axial  flow 
at  Mach  number  M  ,  as  shown  in  Fig.  1.  Any  inflow  distortions  are  viewed  as 
small  perturbations  about  this  undisturbed  state.  The  blades  rotate  with  a 
constant  angular  velocity,  which  in  this  analysis  is  denoted  ,  rather  than 
.a,'  as  used  previously.  In  keeping  with  the  usual  convention,  the  latter 
symbol  will  be  used  for  the  harmonic  time  dependence  below.  Since  the  blade 
boundary  conditions  are  more  easily  expressed  in  blade -fixed  coordinates,  we 
again  express  the  governing  equations  in  these  terms.  In  this  frame,  the 
steady  un disturbed  inflow  has  a  velocity  Un*[U  ♦  i  si  '  J  *,  and  follows 
the  helical  stream  surfaces  defined  by  ^  *  t?  -  j*  *  constant.  The  unsteady 

flow  is  assumed  to  be  a  small  perturbation  about  this  flow,  so  the  disturbance 
field  will  be  irrotation.il  and  isentropic.  The  linearisation  again  allows  ns 
to  apply  the  blade  boundary  conditions  along  the  undisturbed  stream  surfaces, 
so  that  to  first  order  in  the  perturbation  scheme,  blade  thickness  and  camber 
do  not  affect  the  unsteady  loads.  We  also  assume  that  the  undisturbed  relative 
Mach  number  is  subsonic  all  along  the  span. 


1'he  linearised  three-dimensional  equations  of  motion  in  blade-fixed 
coordinates  can  then  be  reduced  to  the  following  single  equation  for  the 


perturbation  pressure,  e  : 
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where  the  operator  3t  is  defined  as  (^~- 
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is  the  force  per  unit  volume  exerted  by  the  Made  row  on  the  fluid,  and  ,1. 
is  the  undisturbed  sound  speed.  (As  will  be  seen  shortly,  the  flow  at  upstream 
Infinity  will  no  longer  be  undisturbed,  due  to  the  presence  of  acoustic  waves. 
Accordingly ,  in  this  section  we  will  use  the  subscript  I  )j  to  denote  condi¬ 
tions  at  the  reference  state.)  The  hard-wall  boundary  condition  at  the  hub  and 
tip,  y>  *  C  ,  can  again  be  expressed  as 
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liquations  t-JOl  and  (U)  are  the  same  as  those  solved  in  the  stead) 
flow  case,  except  for  the  addition  of  the  three  unsteady  terms  at  the  end  of 
(40b)  . 

C.  SOLUTIONS  FOR  PRESSURH  MONOPOLY  AND  PI  POLL  SINGULAR! THIS 

If  we  assume  a  harmonic  c“*’£  time  dependence  for  the  excitation, 
whether  it  be  due  to  blade  flutter  or  inlet  distortion,  the  pressure  and 
velocity  fields  will  also  be  harmonic  at  the  same  frequency  because  of  the 
linearisation  of  the  equations.  Then  the  first  step  in  the  solution  of 
Fqn.  (401  is  to  obtain  the  so-called  pressure  monopole  solution,  .  which 


satisfies : 


jr  6  (.3  -  9t )  6  ( r  -  rt )  <f  <  *  -  e# )  c 1  u'r 


- 


(42) 


subject  to  the  hard-wall  boundary  condition  in  Eqn.  (41). 

This  is  obtained  along  the  same  lines  as  were  followed  for  the  steady  flow 
problem^,  and  the  result  is: 
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where  now, 
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As  before,  3  *  7-M,  a" ;  — ,  and  n  and  k  are  the  azimuthal  and  radial  duct  mode 

'T 

indices,  respectively.  and  are  the  corresponding  radial  eigenfunctions 

and  eigenvalues  which  satisfy  Eqn.  (41). 

In  the  analogous  steady  flow  expression,  Eqn.  (SI)  of  Ref.  1,  it  was 
found  that  since  is  always  greater  than  ^  then  for  subsonic  relative 

Mach  numbers,  (  )“  >  (  )2.  As  a  result,  all  modes  were  below  cut¬ 

off,  with  the  exception  of  the  n  =  o  ,  k  *o  mode.  In  contrast,  a  notable 
feature  of  the  above  solution  is  that  a  given  mode  may  or  may  not  propagate 
undamped  down  the  duct,  depending  on  whether  (  ~%"~  rT  )  is  greater  or  less 
than  (  fi  K  „h  ).  For  modes  below  cut-off,  the  signs  in  the  complex  exponen¬ 
tial  of  Eqn.  (43)  were  chosen  so  that  such  modes  will  decay,  rather  than  grow, 
with  distance  from  the  singularity.  For  propagating  modes,  i.e.,  those  above 
cut-off,  the  choice  is  dictated  by  the  requirement  that  the  group  velocity  of 
the  mode,  which  is  the  velocity  at  which  energy  is  carried  in  the  flow,  always 
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be  away  from  the  singularity.  Alternatively,  a  small  damping  term  can  be  added 
to  the  left  side  of  Eqn.  (42),  which  results  in  even  the  modes  above  cut-off 
having  a  small  real  part  in  the  complex  exponential  of  Eqn.  (43).  The  correct 
sign  choice  is  then  made  by  requiring  that  this  part  always  be  negative,  after 
which  taking  the  limit  of  zero  damping  leads  again  to  Eqn.  (43). 


Right  at  cut-off,  Aot<  =  0  ,  i.e. 


in  xl  +  uJ 


and  Eqn.  (43)  exhibits  singular  behavior.  However,  at  this  condition,  the 
derivation  leading  up  to  (43)  is  no  longer  valid,  so  this  situation  must  be 
treated  as  a  special  case.  The  solution  for  such  a  mode  is, 
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2  -  2.  I  «r>  (4b) 
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where  A  ,  and  3,  are  undetermined  constants.  This  indeterminate  nature  of  the 

"'to 

solution  is  evidently  the  result  of  attempting  to  excite  a  linear  system  at 
one  of  its  resonant  points.  Although  interesting,  such  a  situation  is  not 
likely  to  arise  in  practice.  This  is  because  the  set  of  Knk  in  Eqn.  (45) 
depends  only  on  the  hub/ tip  ratio,  blade  number,  and  inter-blade  phase  angle, 
whereas  the  right-hand  side  is  a  function  also  of  the  rotor  speed,  excitation 
frequency,  and  axial  Mach  number. 


Another  interesting  feature  to  note  is  that  the  n  •  k  *  o  term  in 
Eqn.  (43)  can  be  written  in  the  same  form  as  the  others.  This  was  not  true 
in  the  steady  case,  where  this  mode  had  an  entirely  different  type  of  structure. 
N'or  can  the  correct  steady  term  be  obtained  from  that  in  Eqn.  (43)  by  simply 
setting  to  *  o  .  This  iinomaly  is  explained  by  recognizing  from  Eqn.  (44)  that 
in  the  unsteady  case  \O0  remains  finite  despite  the  fact  that  k  w  *  0  , 

whereas  in  the  steady  problem  \oo  vanishes.  Hence,  in  the  unsteady  case,  this 
mode  is  always  above  cut-off,  but  in  the  steady  case,  it  is  always  right  at 
the  cut-off  condition.  The  correct  form  of  this  mode  in  the  steady  problem 
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should  thus  follow  by  specializing  Eqn.  (46)  to  -u  «  n  «  k  •  0  ,  and  we  have 
shown  that  this  is  indeed  the  case. 

The  next  step  is  to  obtain  an  expression  corresponding  to  a  pressure 
dipole  singularity,  i.e.,  the  solution  to 
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where  r\  is  the  unit  vector  perpendicular  to  the  local  undisturbed  stream  di¬ 
rection,  as  shown  in  Fig.  2.  This  is  obtained  by  differentiating  the  monopole 


solution  along 

this  direction  in 

source  coordinates: 
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Physically,  the  pressure  dipole  solution  represents  the  pressure  field  gener¬ 
ated  by  a  unit  amplitude  harmonic  force  concentrated  at  the  point  (  ra  ,  da  ,  z,  ) 
In  the  next  sub-section,  we  use  it  to  build  up  the  pressure  field  for  the 
entire  blade  row  by  superposition. 


D.  PRESSURE  FIELD  FOR  ENTIRE  BLADE  ROW 

The  solution  to  Eqn.  (40)  for  the  entire  blade  row  is  next  obtained 
by  distributing  the  pressure  dipole  singularities  along  the  undisturbed  stream 
surfaces  ^  in  Fig.  2,  with  an  amplitude  equal  to  the  local  unsteady 

force  exerted  by  the  blade  on  the  flow.  In  general,  we  must  allow  for  a  phase 
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t?rence  between  the  loading  on  successive  bladfsj  denoting  the  local  loading 
on  the  jth  blade  by  A  -©  ,  we  have 
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where  ug  is  the  inter-blade  phase  shift,  and  Af?,  the  loading  on  the  zeroth 
or  reference  blade.  Periodicity  requires  that 
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Higher  values  of  p  would  simply  reproduce  the  same  phase  shifts,  e.g.,  -p  *  S 
is  physically  equivalent  to  -v »  o ,  so  we  need  consider  only  those  values  of  p 
indicated  in  Eqn.  (49b).  Note  however,  that  in  any  given  calculation  -o  , 


and  hence  p.  ,will  assume  a  unique  value. 

ft 


By  superposition,  the  pressure  field  for  the  entire  blade  row  may  be 


written  as 


r  C(r ) 
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where  we  have  set  Apa  =  Ap  .  Substituting  for  -p0  and  ^.8>;ind  transforming  t!ie 
streamwise  integration  to  one  over  the  axial  coordinate  using 

-  f(W  *. 

we  obtain 
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where  now  n  -  +  p 


In  writing  this  result,  it  was  assumed  for  convenience  that  the  axial  projec 

/  AT*  *  * 

tion  of  the  blade  chord,  =  c  /  yi  was  a  constant  independent  of 

radius  (Fig.  2).  Also  note  that  the  same  symbol -p  has  been  used  to  denote 
both  the  unsteady  pressure  and  the  integral  shift  in  acimuthal  mode  number 
caused  by  the  phase  shift  between  blades;  when  used  below  its  meaning  should 
be  clear  fron  the  context  in  which  it  appears. 


UNSTEADY  INTEGRAL  EQUATION 


The  desired  integral  equation  relates  the  unknown  loading,  Je  , 
to  Vn  ,  the  velocity  component  along  the  n  direction  in  Fig.  2.  It  is  this 
velocity  which  is  most  easily  related  to  the  excitation,  whether  inflow  dis 
tortion  or  blade  flutter,  through  the  surface  tangency. condition .  The  ex¬ 
pression  for  v„  may  be  derived  from  the  normal  component  of  the  momentum 
equation : 


which,  for  harmonic  time  dependence  reduces  to 


where  it  is  understood  that  the  variables  now  refer  only  to  their  spatial 
amplitudes.  This  equation  may  be  integrated  from  upstream  infinity  along  the 
undisturbed  stream  direction,  S  •  The  constant  of  integration  which  must  be 
supplied  is  the  value  of  ir0  at  upstream  infinity.  Following  the  same  line  of 
reasoning  given  after  Eqn.  (44),  we  imagine  that  a  small  amount  of  dissipation 
is  present  in  the  system  while  performing  the  integration,  which  insures  that 
lr„  (-»  )  *  O  .  The  limit  of  :ero  dissipation  is  then  taken  in  the  result, 
which  reads 
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where  the  integration  is  to  be  tarried  out  with  r  and  „  (or  ?  )  held  constant 

The  normal  derivative  of  the  pressure  is  obtained  from  Eqn.  (51)  as 
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where  d\x>  and  Hex 


respectively,  and 


'are  the  Dirac  delta  function  and  Heaviside  step  function. 


f  (d  ~  —  _  /m_y  /  nn  ^ 

r  (fiH  u)  (~u~)  (55, 

The  force  F8  is  considered  positive  when  it  points  in  the  +n 
direction;  for  the  entire  blade  row  it  can  be  expressed  as. 
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Whom  the  minus  sign  in  front  is  required  by  the  fact  that  in  our  convention 
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Using  this  and  the  definition  of  ixd  in  Eqn.  (49b),  we  obtain 


Ap  C'*.  4 


Substituting  (54)  and  (56)  into  (53)  gives 
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where  we've  implicitly  assumed  0  &  2  $  ,  so  that  the  Heaviside  step  func¬ 
tions  in  the  first  term  could  be  set  to  unity.  Since  r,  (, ,  zt  and  £  are  eac 
constant  during  the  integration  over  £' ,  the  latter  can  be  done  analytically 
After  some  manipulation,  the  result  is: 
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where  f(r>  was  defined  in  (55)  and 

n  /iif  \  /  niltw  v  ,co. 
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Consider  the  convergence  of  the  rv>  and  k  series  in  Eqn.  (58).  In 
the  last  two  terms,  convergence  of  the  series  is  assured  by  the  presence  of 
the  factor  e  oK  i  *‘l  .  This  is  because  for  m  (and  hence  n  )  fixed,  as 
k  becomes  large,  K„k  grows  linearly  with  k  A 

~  77a  k  >  >  n  ,  n  fixed  (bO) 

and  beyond  a  certain  point,  the  modes  cut-off  (Eqn.  (44)).  From  then  on,  A  k 

remains  real  and  grows  with  k  ,  thus  insuring  exponential  convergence. 

Conversely,  in  the  limit  of  large  )  m)  (and  hence  |o)  ),  for  k  fixed, 
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it  is  known  that  "  : 
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Since  we  have  assumed  /3  ‘  >  0  ,  convergence  is  again  assured. 

Unfortunately,  this  exponential  decay  is  not  operative  in  the  first 

group  of  terms  under  the  double  summation,  i.e.,  those  proportional  to  H 

these  diverge  with  respect  to  both  m  and  k.  The  first  tenn  in  Eqn.  (SS)  also 

poses  a  problem,  in  that  it  exhibits  singular  behavior  on  the  blade  surfaces 
2  r 1 1 

Zj  =  ~af~~  >  which  is  where  we  ultimately  wish  to  evaluate  it.  Similar 
difficulties  arose  in  the  steady  flow  problem;  to  remove  them  in  the  present 
case,  a  procedure  very  much  analogous  to  that  used  in  Ref.  1  was  followed. 
Namely,  the  same  group  of  divergent  terms  was  added  and  subtracted  from  the 
original  divergent  series.  The  added  terms  can  ultimately  be  summed  in  the 
context  of  generalized  functions  to  give  a  sequence  of  delta  functions  which 
exactly  cancels  the  first  term  of  Eqn .  (58).  The  original  minus  the  sub¬ 
tracted  terms,  through  use  of  the  fact  that  each  of  the  satisfies 

Bessel's  equation,  can  be  rewritten  in  terms  of  the  first  and  second  radial 
derivatives  of  the  R„u  •  If  two  successive  integrations  by  parts  in  the 
radial  coordinate  are  then  made,  and  advantage  taken  of  the  fact  that  both 
Ap  and  R^ satisfy  Eqn.  (41),  the  result  is  a  series  which  is  convergent  with 
respect  to  both  m  and  k. 


This  group  of  terms,  as  in  the  steady  flow  case1,  represents  the 
wakes  shed  by  the  blades.  In  the  present  case,  the  wakes  contain  both  trail¬ 
ing  vorticity  produced  by  radial  gradients  in  the  circulation,  and  shed  vor- 
ticity  created  by  the  fact  that  the  circulation  is  now  varying  in  time  as  well 
Since  the  manipulations  described  above  are  quite  similar  to  those  used  in 
Ref.  1,  but  somewhat  lengthier,  they  are  not  reproduced  here.  The  end  result 
for  v„  is: 
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where  it  should  be  recalled  that  n  r>-.i5»©and  \  ,  ,  f ^rj  and.7<r)are  defined 
in  (44),  (55)  and  (59),  respectively. 

The  above  expression  contains  only  convergent  series.  The  first 
two  lines  are  the  contributions  to  the  wake  from  the  axisymmetric  modes,  n  = 

0.  From  the  relation  between  n  and  we  see  that  such  modes  can  arise  only 
from  the  w  =  0  term,  and  then  only  in  the  special  case  -p  =  0  i.e.,  ;ero  phase 

lag  between  blades.  Accordingly,  these  are  written  with  the  Kronecker  delta 
symbol,  6  p<>  ,  which  is  one  when  p  =  0  and  zero  otherwise.  The  next  eight 
lines  are  the  contributions  to  the  wake  terms  from  the  |nl  >  0  modes.  All 
of  the  wake  terms  contain  the  exponential  factor  e  ^  u(z  -*)  ^ 

as  well  as  the  implied  C  time  dependence.  Along  the  undisturbed 

helical  streamlines  q  =  constant),  this  can  be  rewritten  as 

[t-  ‘^4s2] 


Thus,  an  observer  moving  in  this  direction  with  velocity  IL  senses  a  frozen 
velocity  pattern  from  these  terms.  This  reflects  the  fact  that,  in  a  linear- 
iced  analysis,  the  wakes  convect  downstream  unattenuated  at  the  freest  ream 
relative  velocity. 


The  last  four  lines  in  Eqn.  (65)  represent  the  wave  field  that  pro¬ 
pagates  away  from  the  rotor  in  both  directions.  This  consists  of  both  damped 
and  undamped  waves,  or  modes;  whether  a  particular  mode  is  damped  or  not  is 
determined  by  the  cut-off  criterion  described  by  Eqn.  (44).  The  modes  above 
cut-off  represent  a  loss  of  energy  through  acoustic  radiation.  Hence,  a  useful 
by-product  of  the  linearized  analysis  is  that  it  will  also  predict  the  total 
flux  of  acoustic  energy,  as  well  as  its  distribution  among  the  various  propa¬ 
gating  duct  modes. 


The  only  steps  remaining  in  the  formulation  of  the  unsteady  integral 
equation  are  to  specialize  Eqn.  (65)  to  the  reference  blade  surface,  r,»  o, 
and  to  specify  the  normal  velocity  there,  ( r,  t,  *c7,  l  ),  in  terms  of  some 
prescribed  excitation.  The  latter  is  accomplished  by  requiring  that  at  each 
instant  the  streamlines  remain  tangent  to  the  blade  surface.  The  excitation 
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could  be  an  inflow  distortion  pattern,  such  as  might  result  from  the  wakes  of 
an  upstream  blade  row  or  distortions  in  the  engine  inlet  itself.  Or,  v„  could 
be  considered  to  be  the  result  of  vibrations  of  the  blades  themselves,  i.e., 
the  so-valled  flutter  problem. 

In  any  event,  once  the  left-hand  side  of  Eqn.  (65)  is  prescribed, 
the  integral  equation  which  results  must  be  inverted  to  obtain  the  unknown 
loading  amplitude,  A-p  .  Generally  speaking,  at  each  point  there  will  be  a 
phase  shift  between  the  unsteady  loading  and  the  normal  velocity  which  will 
depend  on  position,  reduced  frequency,  and  Mach  number.  This  was  not  true  in 
the  steady  flow  problem,  where  by  definition,  everything  is  "in  phase."  To 
account  for  this.Ap  must  now  be  allowed  to  be  a  complex  quantity,  having  both 
real  and  imaginary  parts.  What  we  actually  have  to  solve  then,  is  two  integral 
equations,  one  each  from  equating  the  real  and  imaginary  parts  of  Eqn.  (63). 

Equation  (63)  is  far  too  complicated  to  hope  for  an  analytical  solu¬ 
tion,  and  so  a  numerical  inversion  procedure  must  be  developed.  At  present, 
it  appears  an  extension  of  the  procedure  applied  in  Section  I  to  the  steady 
problem  should  work  here  as  well.  Specifically,  Ap  would  be  expanded  in  the 
same  finite  series  as  in  Eqn.  (12),  except  that  now  the  unknown  coefficients 
would  be  complex.  In  addition  to  automatically  satisfying  the  appropriate 
leading  and  trailing  edge  conditions,  use  of  this  expansion  will  again  lead  to 
greater  efficiency  in  computing  the  axial  and  radial  integrations  in  (63)  .  The 
axial  integrations  will  be  of  the  same  form  as  in  Eqn.  (14),  although  the  form 
of  the  constant  parameters  will  change  somewhat.  Thus,  the  analytical  evalua¬ 
tion  procedure  given  in  Appendix  C  can  again  be  used  to  good  advantage. 

The  majority  of  terms  in  Eqn.  (63)  will  also  lead  to  radial  inte¬ 
grations  of  the  same  form  as  in  (13).  The  exceptions  are  those  portions  of  the 
wake  terms  which  involve  reciprocal  powers  of  [  i  +  ];  however,  such 

factors  vary  monotonical ly  from  hub  to  tip,  and  so  should  be  expressible  as  a 
finite  series  of  integral  powers  of  <T0  =  with  a  high  degree  of  accuracy. 

Hence,  we  believe  that  the  radial  integrations  could  again  be  performed  an¬ 
alytically,  as  outlined  in  Appendix  B. 

I 

« 


46 


If  NX  x  NJ  terms  are  retained  in  the  expansion  for  So  ,  we  would 
have  a  total  of  2  x  NI  x  NJ  scalar  unknowns.  Substituting  this  expansion  into 
the  integral  equations  described  above,  and  evaluating  the  resulting  set  of 
equations  at  an  equal  number  of  collocation  points,  will  yield  a  system  of 
2  x  NI  x  NJ  simultaneous,  linear,  algebraic  equations.  These  can  again  be 
viewed  as  a  single  matrix  equation,  with  the  elements  of  the  coefficient  matrix 
being  expressed  in  terms  of  the  axial  and  radial  integrations  discussed  above. 
The  inversion  of  the  matrix  equation  could  then  be  carried  out  in  a  straight¬ 
forward  manner  using  standard  techniques. 


F.  CONCLUSIONS 

A  linearited  three-dimensional  analysis  has  been  developed  for 
studying  the  unsteady  loading  on  an  annular  blade  row  in  compressible  flow. 
Pressure  dipole  singularities  are  used  to  represent  the  blades  in  a  lifting- 
surface  type  analysis.  The  result  is  an  integral  equation  which  relates  the 
unknown  loading  to  the  prescribed  variations  in  normal  velocity  at  the  blade 
surfaces.  The  latter  can  be  the  result  of  either  an  inflow  distortion  pattern 
or  vibrations  of  the  blades  themselves.  Terms  can  be  identified  in  the  inte¬ 
gral  equation  which  represent  the  blade  wakes.  These  include  elements  of  both 
trailing  vorticity,  produced  by  spanwise  loading  gradients,  and  shed  vorticity, 
which  results  from  the  temporal  fluctuations  in  loading.  The  wakes  can  be 
shown  to  convect  downstream  along  the  helical  undisturbed  stream  surfaces. 

Also  included  are  terms  which  represent  waves  propagating  unattenuated  away 
from  the  blade  row  in  both  directions.  This  acoustic  radiation  field  would  be 
predicted  as  a  by-product  of  the  solution,  with  regard  to  both  total  acoustic 
energy  flux,  and  its  distribution  among  the  various  propagating  duct  modes. 

The  unsteady  integral  equation,  though  somewhat  more  complicated 
than  the  analogous  steady  flow  equation  discussed  in  Section  I,  contains  axial 
and  radial  integrations  of  the  same  basic  analytical  form.  This  suggests  that 
a  suitable  extension  of  the  numerical  inversion  procedure  developed  for  the 
steady  case  should  work  here  as  well,  and  this  is  the  approach  that  is  being 
pursued . 


SECTION  III 


FURTHER  COMPARISONS  OF  APPROXIMATE  ROTOR-STATOR  INTERACTION  THEORY 
WITH  ACOUSTIC  DATA  FROM  ANNULAR  CASCADE 

A.  INTRODUCTION 

Under  the  previous  contract,  an  approximate  theoretical  model  was 
developed  to  predict  the  noise  generated  by  rotor-stator  interaction.'  In  order 
to  assess  the  accuracy  ot'  this  theoretical  model,  experimental  acoustic  data 
were  obtained  in  an  annular  cascade  facility  upstream  of  a  rotor- stator  pair, 
and  compared  with  the  theoretical  predictions.  Unfortunately,  power  constraints 
on  the  rig  restricted  the  acoustic  tests  to  low  rpm  where  only  the  fourth  and 
higher  harmonics  could  propagate  undamped  in  the  duct.  The  prediction  of  these 
tones  was  felt  to  be  a  rather  severe  test  of  the  theory.  While  the  theory  was 
able  to  predict  cut-off  frequencies  and  follow  relative  trends  in  the  data 
adequately,  it  significantly  underestimated  absolute  levels. 

Additional  numerical  studies  using  the  theory  showed  that  the  prediction 
of  such  high  harmonics  is  critically  dependent  on  the  assumed  shape  for  the  rotor 
blade  wakes.  Very  small  changes  in  assumed  wake  shape  can  lead  to  large  changes 
in  the  predicted  noise  levels,  particularly  for  the  higher  harmonics.  The 
approximate  theory  relies  on  isolated  airfoil  data  for  the  wake  shapes.  It  may 
be  that  the  prediction  of  interference  noise  at  these  high  harmonics  of  blade 
passage  frequency  cannot  be  accomplished  in  any  practical  sense  until  we  are 
able  to  predict  the  wake  structure  to  a  much  higher  degree  of  accuracy. 

On  the  other  hand,  because  of  the  much  higher  rpm  at  which  typical 
aircraft  fans  and  compressors  operate,  it  is  the  fundamental  blade  passage 
frequency,  and  perhaps  its  second  harmonic,  which  account  for  most  of  the  noise 
problem.  At  such  speeds,  the  higher  harmonics  fall  outside  the  range  of  human 
hearing.  The  theory  indicated  that  the  first  few  harmonics  of  blade  passage 
frequency  are  relatively  insensitive  to  the  shape  of  the  rotor  blade  wakes  and 
that  the  level  of  the  interaction  noise  at  these  harmonics,  if  excited,  will  be 
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much  larger  than  the  levels  of  the  higher  harmonics.  Therefore,  it  was  felt 
that  a  more  reasonable  test  of  the  theory  would  be  to  compare  its  predictions 
against  data  taken  for  the  iower  harmonics.  This  required  modifying  the  drive 
mechanism  in  the  annular  cascade  so  that  higher  speeds  could  be  achieved.  The 
repowered  rig  and  the  experimental  procedures  are  described  in  Section  III-B, 
followed  by  comparisons  with  the  theoretical  predictions  in  Section  III-C. 

B.  ROTOR-STATOR  INTERACTION  NOISE  MEASUREMENTS 

1 .  Rotating  Annular  Cascade  Facility 

As  part  of  the  work  under  a  previous  program,  Contract  AF35  (615) -335 
an  annular  cascade  facility  was  designed  and  fabricated.  Its  principal  purpose 
is  to  provide  fundamental  experimental  data  during  and  prior  to  the  occurrence 
of  rotating  stall  in  order  to  improve  our  understanding  of  the  phenomena  and 
for  use  as  a  guide  in  improving  the  theoretical  rotating  stall  analysis.  In 
addition  to  the  study  of  rotating  stall,  the  facility  has  also  been  used  to 
provide  acoustic  data  for  comparison  with  theory.  A  detailed  description  of  the 
annular  cascade  facility  has  been  presented  in  Ref.  31.  Only  a  brief  descrip¬ 
tion  of  the  facility  is  given  here. 

The  annular  cascade  facility  consists  of  a  test  section  built 
around  the  outer  front  casing  of  a  J-79  jet-engine  compressor  with  a  Calspan 
fabricated  hub.  The  facility  includes  a  bell-mouth  inlet  on  the  outer  casing  and 
a  bullet  nose  on  the  hub  to  provide  a  smooth  flow  of  air  to  the  test  section. 
Outlet  ducting  is  connected  to  an  independently  variable  source  of  suction  to 
provide  the  required  flow  through  the  annulus.  An  electrically  powered  two- 
speed  axial  flow  fan  is  used  as  the  source  of  suction.  Continuous  control  of 
the  mass  flow  is  achieved  through  the  use  of  variable  inlet  guide  vanes  to  the 
fan  and  a  variable  damper  in  the  fan  exit  flow. 


The  configuration  of  the  annular  cascade  used  for  the  rotor-stator 
interaction  studies  is  shown  in  Fig.  14.  The  test  section  forms  a  circular 
annulus  with  an  outer  diameter  of  29.35  inches  and  an  inner  diameter  of  25.55 
inches  which  provides  a  hub-to-tip  ratio  of  0.S0.  The  rotor-stator  stage 
studied  is  designated  Rotor  Set  No.  1  (.46  blades)  and  Stator  Set  No.  1  i54 
vanes)  in  Ref.  31.  These  are  modified  blade  rows  from  the  fifth  stage  of  a 
J-~9  compressor.  Their  characteristics  and  performance  when  used  as  isolated 
blade  rows  in  the  annular  cascade  have  been  presented  in  Ref.  51.  The  semi¬ 
chords  of  the  rotor  and  stator  were  0.0604  ft.  and  0.054  ft.,  respectively,  and 
the  axial  separation  of  the  mid-chord  planes  was  0.125  ft.  The  mid-annulus 
stagger  angle  of  the  rotor  was  40  degrees  in  all  tests.  Two  different  stator 
stagger  angles  (. 3 ~ .  2  and  2S.2  degrees)  were  tested. 

The  rotor  in  the  annular  cascade  is  driven  by  a  hydraulic  motor.  An 
external  hydraulic  pump  system  is  used  to  provide  power  for  the  hydraulic 
motor.  In  the  tests  of  Ref.  1,  the  hydraulic  pump  system  was  powered  by  a 
30  horsepower  electric  motor.  Kith  this  system,  the  maximum  rotor  speed 
attainable  was  1450  rpm.  This  speed  was  insufficient  to  allow  propagation 
of  sound  at  the  fundamental  blade  passage  frequency.  It  was  calculated  that  an 
increase  in  rotor  speed  from  the  maximum  of  1450  rpm  attained  in  the  past 
experiments  to  approximately  1700  rpm  should  be  sufficient  to  exceed  the  cut¬ 
off  condition  for  propagation  of  the  fundamental  blade  passage  tone  in  the  rotor 
stator  stage.  Since  the  time  that  the  past  experiments  were  performed,  the 
rotating  annular  cascade  was  modified  to  include  a  50  horsepower  hydraulic  pump 
system  in  the  rotor  drive,  and  changes  were  made  to  reduce  pressure  losses  in 
the  hydraulic  valve  system  used  for  rotor  speed  control.  With  these  changes, 
a  maximum  rotor  speed  of  1"60  rpm  was  attained.  Acoustic  tests  in  the  modified 
annular  cascade  succeeded  in  exciting  the  fundamental  as  well  as  higher  harmonic 
of  blade  passage  frequency. 

2.  Acoustic  Measurements 


The  Sound  Pressure  Levels  produced  by  rotor-stator  interaction 
were  measured  on  the  outer  casing  upstream  of  the  rotor.  In  taking  these  data, 
the  fan  system  downstream  of  the  annular  cascade  was  turned  off  and  the  fan  was 


allowed  to  rotate  freely  under  the  influence  of  the  flow  generated  by  the  rotor 
in  the  annular  cascade.  In  addition,  all  dampers  in  the  downstream  drive  system 
were  opened  wide.  It  was  convenient  to  use  this  wide  open  configuration  because 
it  was  found  to  provide  nearly  constant  inlet  swirl  angles  relative  to  the  rotor. 
That  is,  the  mean  axial  velocity  in  the  annular  cascade  increased  in  direct 
proportion  to  the  rotor  angular  velocity.  This  feature  made  application  of  the 
rotor-stator  interaction  theory  much  simpler  because  the  dimensionless  steady 
state  loss  and  turning  performance  data  for  the  rotor  also  remained  approximately 
constant  with  changing  rotor  speeds. 

The  instrumentation  used  for  the  acoustic  tests  is  shown  schemati¬ 
cally  in  Fig.  15.  A  Bruel  and  Kjaer  (BSK)  condenser  microphone  (Type  4135) 
mounted  with  its  diaphragm  flush  with  the  inner  surface  of  the  compressor  casing 
(Fig.  14)  was  used  to  measure  the  rotor-stator  interaction  noise.  The  micro¬ 
phone  was  powered  by  a  BSK  Type  2604  microphone  amplifier.  The  output  from  the 
amplifier  was  processed  on  line  by  a  General  Radio  Type  1921  real  time  analyzer 
to  obtain  one-third  octave  spectra  of  the  signals.  Absolute  calibration  of  the 
complete  system  was  obtained  by  recording  a  114  db  signal  from  a  General  Radio 
Type  1562-A  sound  level  calibrator  on  the  recorder  chart  of  the  spectrum  analyzer. 
The  calibration  signal  was  recorded  at  the  start  of  a  test  series  and  again  at 
the  end  of  a  test  series  to  check  for  drift  of  the  system.  The  observed  drift 
was  0.5  db  or  less  in  all  tests. 

The  summed  output  from  the  multifilter  in  the  GR  real  time  analyzer 
was  displayed  on  a  dual  beam  oscilloscope  along  with  a  reference  signal  from  an 
RC  oscillator.  The  multifilter  summed  output  was  conditioned  by  attentuating 
all  one  third-octave  levels  with  center  frequencies  of  500  Hz  and  below.  This 
was  done  to  reduce  noise  at  frequencies  less  than  the  fundamental  blade  passage 
frequency.  Oscilloscope  records  of  the  conditioned  noise  signal  are  shown  in 
Fig.  16.  In  these  photographs,  the  lower  trace  is  the  noise  signal  and  the 
upper  trace  is  the  output  from  the  RC  oscillator  tuned  to  the  same  frequency 
as  the  predominant  noise  harmonic.  The  jitter  in  the  noise  signals  is  caused 
by  the  presence  of  frequencies  other  than  the  predominant  harmonic.  Tuning  of 
the  audio  oscillator  was  accomplished  by  setting  the  oscilloscope  trigger 
mechanism  to  fire  on  the  audio  oscillator  signal  and  then  adjusting  the  oscillator 


frequency  until  a  stationary  noise  signal  was  attained.  In  some  tests,  the 
correct  frequency  was  difficult  to  determine  because  of  the  iitter  in  the  noise 
signal.  The  correct  frequency  was  then  determined  by  forming  a  Lissaious 
figure  from  the  noise  and  oscillator  signals,  and  tuning  the  oscillator  frequency 
until  a  stationary  Lissaious  figure  was  obtained.  The  frequency  was  then  read 
on  a  Berkley  universal  counter  and  timer.  A  similar  counter  was  used  to  deter¬ 
mine  the  rotor  speed  from  signals  generated  by  a  magnetic  pickup  and  120  tooth 
gear  on  the  rotor  drive  shaft. 

The  procedure  described  in  the  previous  paragraph  was  used  to  determine 
the  frequency  of  the  predominant  pure  tone  noise  for  all  tests  where  a  strong 
blade  row  Interaction  signal  was  obtained.  Thus,  it  was  possible  to  identify 
accurately  which  harmonic  of  rotor  blade  passage  frequency  was  being  excited 
by  the  interaction  mechanism.  In  the  upper  photograph  of  Fig.  16,  the  predomi¬ 
nant  frequency  is  the  fundamental  of  blade  passage  frequency.  In  the  lower 
photograph,  both  the  fundamental  and  second  harmonic  of  blade  passage  frequency 
are  apparent . 

Samples  of  the  one-third  octave  spectra  of  the  interaction  noise  are 
shown  in  Figs.  1"  and  IS.  These  figures  are  direct  reproductions  of  the  output 
from  the  spectrum  analyser.  Absolute  Sound  Pressure  Level  for  any  one-third 
octave  band  is  obtained  by  correcting  the  relative  db  levels  for  background 
noise,  system  gain,  microphone  response  and  filter  characteristics.  The  system 
gain  correction  is  independent  of  frequency  and  is  noted  on  the  lower  left  hand 
side  of  each  figure.  For  the  particular  example  of  Fig.  1',  system  gain  is 
corrected  for  by  adding  80  db  to  all  relative  levels.  For  Fig.  IS.  the  gain 
correction  is  80.5  db.  The  remaining  corrections  will  be  discussed  shortly. 

Figures  1'a  through  l'c  have  been  chosen  to  illustrate  the  extremely 
sharp  rise  in  interaction  noise  as  rotor  speed  is  increased  over  a  very  small 
range  of  rpm  from  below  cutoff  to  above  cutoff  of  the  fourth  rotor  blade 
passage  harmonic.  Figure  l"a  corresponds  to  conditions  slightly  below  cutoff. 
Here  the  noise  spectrum  is  generated  primarily  by  turbulent  pressure  fluctuations 
within  the  boundary  layer  on  the  outer  wall  and  by  the  free-wheeling  fan  down¬ 
stream  of  the  annular  cascade  test  section.  Figs,  l'h  and  l'c  show  the  sharp 
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rise  in  fourth  harmonic  interaction  noise  as  the  rotor  rpm  increased  to  above 
cutoff  conditions.  This  rise  can  be  seen  in  the  third  octave  bands  centered 
at  5150  and  4000  He .  Both  of  these  third  octave  bands  respond  because  the  pure 
tone  frequency  is  nearly  mid-way  between  these  two  bands.  The  calculated  fre¬ 
quency  of  the  fourth  harmonic  is  indicated  by  an  arrow  on  the  spectra;  the  numer¬ 
ical  value  is  shown  on  the  lower  right  hand  side  of  Figs,  l'a  through  l“c. 

The  corresponding  measured  pure  tone  frequency  is  also  indicated  on  Figs,  l'b 
and  l"c. 

Figs.  ISa  through  lSd  are  spectra  obtained  at  rotor  speeds  sufficient 
to  excite  the  fundamental  blade  passage  frequency.  Here  the  fundamental  fre¬ 
quency  and  the  second  and  fourth  harmonics  are  apparent  in  the  spectra.  In 
Fig.  18a  the  fundamental  frequency  is  predominant  and  smaller  Sound  Pressure 
Levels  are  evident  at  the  second  and  fourth  harmonic  frequencies.  When  the  rotor 
speed  is  increased  by  five  rpm  iFig.  iSb) ,  the  second  harmonic  disappears  into 
the  background  noise.  When  the  rotor  speed  is  again  increased  by  small  amounts 
Figs.  15c  and  lSdt,  the  second  harmonic  grows  rapidly  while  the  fundamental 
harmonic  decays.  The  oscilloscope  records  in  Fig.  lo  correspond  to  the  spectra 
presented  in  Figs.  ISa  and  18c.  The  rapid  growth  and  decay  of  Sound  Pressure 
Levels  at  various  harmonics  of  blade  passage  frequency  with  small  changes  in 
rotor  speed  was  observed  throughout  the  complete  test  series. 

As  noted  earlier,  the  relative  Sound  Pressure  Levels  obtained  from 
spectra  such  as  those  shown  in  Figs.  1'  and  18  require  several  corrections  in 
addition  to  the  correction  for  system  gain.  These  are  as  follows: 

(a)  Background  Level  Correction  -  The  Sound  Pressure  Levels  in  the 

third-octave  bands  contain  contributions  from  the  background  noise 
levels  as  well  as  the  pure  tone  noise  from  rotor-stator  inter¬ 
action.  The  background  noise  level  was  estimated  by  inspection  of 
the  third-octave  bands  adjacent  to  those  responding  to  the  various 
harmonics  of  blade  passage  frequency.  An  example  of  the  procedure 
is  shown  in  Fig.  18a.  In  this  case,  the  relative  background 
level  for  the  third-octave  band  centered  at  1250  He  is  estimated 


by  drawing  a  straight  line  between  the  spectrum  levels  centered  at 
1000  Hz  and  2000  Hz.  This  line  intersects  the  1250  Hz  band  at  a 
spectrum  level  of  51.5  db.  The  spectrum  level  at  the  third-octave 
band  centered  at  1600  Hz  was  not  used  in  this  case  because  it 
contains  some  contribution  from  the  fundamental  blade  passage 
frequency.  Ignoring  the  spectrum  level  at  1600  Hz  makes  very 
little  difference  to  the  background  correction  in  this  case.  If 
the  measured  spectrum  level  is  designated  5^  and  the  background 
level  8m  ,  then  the  corrected  spectrum  level  is  given  by 


SPL  Corrected  for  Background  = 
or  equivalent lv 

Background  Correction  = 


rc>X Oq  i  l^a  ' 

o  L  f  70 


70  J 


■(Adb  /  7  o) 


where 


Adb  »  5m  -  Bm  db 


The  background  correction  is  shown  as  a  function  of  Adb  in  Fig. 
19,  which  is  the  same  as  Fig.  6-5  of  Ref.  52. 


(b)  Filter  Correction  -  At  constant  rotor  speed,  the  rotor-stator 
interference  noise  at  the  various  harmonics  of  blade  passage 
frequency  are  essentially  pure  tones  at  specific  frequencies. 

If  the  frequency  of  these  pure  tones  does  not  correspond  to  one 
of  the  center  frequencies  in  the  third-octave  filters  of  the 
real  time  analyzer,  then  the  measured  third-octave  level  will  be 
lower  than  the  pure  tone  level  by  an  amount  depending  on  the 
ratio  between  the  pure  tone  frequency  and  the  third-octave  center 
frequency.  This  arises  because  of  the  characteristic  shape  of  the 
third-octave  filter  curve.  The  measured  response  of  the  filters 
in  the  real-time  analyzer  used  in  these  tests  is  shown  in  Fig.  20. 
This  curve  is  an  average  of  the  results  obtained  for  filters  with 
three  different  center  frequencies,  .  The  measured  noise 
data  for  a  harmonic  frequency,  f  t  were  corrected  by  using  this 
curve  to  determine  the  attentuation  of  the  filter  at  the  appropriate 
frequency  ratio,  f  /f  .  IVhen  two  third-octave  bands  responded  at 
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approximate  1  >•  equal  levels,  such  as  shown  m  l-'  i  >j  s  .  l~b  and  l  ~c  , 
both  values  of  measured  third-octave  levels  were  used  sepuratelv 
to  estimate  the  correct  pure  tone  sound  pressure  level  uisinr. 
the  correct  values  of  f  f  for  each  hand).  In  these  cases,  the 
two  answers  generally  agreed  to  within  1  db  or  better. 

(c)  Microphone  Response  Correction  -  The  frequenc)  response  of  t he 

B;',k  condenser  microphone  was  not  completely  flat  in  the  frequency 
range  above  1000  Ht.  This  is  shown  in  Fig.  Jl.  The  lower  curve 
in  this  figure  is  the  appropriate  curve  for  the  microphone  flush- 
mounted  in  the  annular  cascade  without  its  protection  grid.  The 
overall  system  gain  was  set  by  using  the  111  db  i!R  sound  level 
calibrator  at  a  frequency  of  1000  H;.  Thus  the  lower  curve  in 
Fig.  -1  can  be  used  directly  to  correct  the  results  for  micro¬ 
phone  response  at  the  harmonic  frequencies  of  interest. 

Numerical  examples  of  the  data  reduction  technique  are  given  in  Tables 
.1  and  a  for  the  spectra  presented  in  Figs.  ISa  and  l"c.  The  example  in  Table  f 
is  for  a  spectrum  in  which  several  harmonics  of  blade  passage  frequency  are 
present  in  separate  third-octave  bands.  The  example  in  Table  a  is  for  a  spectrum 
in  which  two  adjacent  third-octave  bands  are  responding  at  nearly  equal  levels 
to  the  same  harmonic  of  blade  passage  frequency.  In  this  case,  the  data  for 
both  third-octave  bands  have  been  analysed  to  show  the  possible  error  when 
the  pure  tone  frequency  is  not  near  any  of  the  center  frequencies.  The  results 
from  separate  analysis  of  the  data  from  the  two  third-octave  filters  are  1JS.1  db 
and  lJ'.a  db.  This  is  within  the  1  db  difference  mentioned  previously  for 
such  cases. 

The  results  of  the  interaction  noise  measurements  are  presented  in 
Tigs.  22  through  JP.  inch  figure  shows  the  corrected  STl.  of  a  particular  blade 
passage  harmonic  as  a  function  of  rotor  speed.  The  approx im.ute  background  noise 
levels  are  also  shown  in  each  figure.  The  small  downward  facing  arrows  in  these 
figures  indicate  points  where  the  spectrum  levels  were  low  enough  to  be 
obscured  by  the  background  level.  Spectrum  levels  within  o  db  of  the  background 
contain  background  level  corrections  which  exceed  1  db  and  should  be  considered 
as  higlilv  approximate  in  view  of  the  method  used  to  estimate  the  background  noise 


SAMPLE  OF  ACOUSTIC  DATA  REDUCTION  FOR 
SPECTRUM  WITH  MULTIPLE  HARMONICS  OF  BLADE  PASSAGE  FREQUENCY 


Rotor  RPM  =  1695. S  (See  Fig.  13a) 


1695.5 

Blade  Passage  Frequency  (46  blades)  =  46  x  60 


1500  He 


Blade  Passage  Harmonic 

1 

2 

5 

Harmonic  Frequency,  f(Hz) 

1500 

2600 

5200 

Third-Octave  Center  Frequency,  f  (Hz) 

1250 

2500 

5000 

£ 

Frequency  Ratio,  /fo 

1.040 

1.040 

1.040 

Measured  Spectrum  Level,  SN)  (Relative  db) 

52.2 

55.0 

54.2 

Estimated  Background  Level,  BV|  (Relative  db) 

51.5 

50.7 

29.7 

-  BV|  (db) 

20.7 

4 . 5 

4.5 

Background  Correction  (db)  (Fig.  19) 

0 

-2.0 

-1.9 

S\|  Corrected  for  Background  (Relative  db) 

52.2 

55 . 0 

52.5 

System  Gain  Correction  (db) 

+80.5 

+80.5 

+80.5 

Filter  Correction  (db)  (Fig.  20  for  ^'/f  =1.040) 

+  0.2 

+  0.2 

+  0.2 

Microphone  Response  Correction  (db)  (Fig.  21) 

+  0.2 

+  0.6 

+  1.S 

Final  Corrected  Sound  Pressure  Level  (db) 

155.1 

114.5 

114.8 

Table  3 

SAMPLE  OF  ACOUSTIC  DATA  REDUCTION  FOR  SPECTRUM  WITH 
TWO  ADJACENT  THIRD-OCTAVE  BANDS  RESPONDING  NEARLY  EQUALLY  TO 
THE  SAME  HARMONIC  OF  BLADE  PASSAGE  FREQUENCY 

Rotor  RPM  =  1160.5  (See  Fig.  l"c) 

1160.5 

Blade  Passage  Frequency  L 4c>  blades)  =  46  x  60  =  889."  H: 

Fourth  Harmonic  of  BPF,  f  =  3559  Hz 


Third  Octave  Center  Frequency,  f  (Hz ) 

3150 

4000 

Frequency  Ratio,  / 1' 

0 

1 .  1 30 

0 . 890 

Measured  Spectrum  Level,  S„,  (Relative  db) 

41.6 

43.0 

Estimated  Background  Level,  BV),  (Relative  db) 

25." 

25.  5 

^M  "  BM 

17.9 

19.  " 

Background  Correction  (db)  Fig.  19) 

-0.  1 

0 

S  Corrected  for  Background  (Relative  db) 

41.5 

43.0 

System  Gain  Correction  (db 

+  80.0 

+  80.0 

Filter  Correction  (db)  (Fig.  20) 

+  5.5 

+  5.2 

Microphone  Response  Correction  (db)  (Fig.  211 

+  1.  1 

+  1.  1 

Final  Corrected  Sound  Pressure  Level  (db) 

128.  1 

12".  5 

57 


Vll  ot'  the  acoustic  measurements  were  made  with  t fie  rotor  stagger  angle 
held  constant  at  40  degrees.  However,  two  stator  stagger  angles  were  used. 

Figs.  2d  through  25  show  the  data  measured  with  a  stator  stagger  angle  of  5". 2 
degrees  while  Figs.  2t>  through  29  show  similar  data  with  a  stator  stagger  angle 
of  2S.2  degrees.  Data  were  obtained  for  the  fundamental  blade  passage  frequency 
and  the  second,  fourth  and  fifth  harmonics.  The  third  harmonic  was  detected  in 
only  two  of  the  195  spectra  obtained  in  the  complete  test  series.  These  two 
data  points  are  shown  in  Fig.  2"  along  with  the  second  harmonic  da’a  measured 
with  a  stator  stagger  angle  of  28.2  degrees.  There  is  evidence  that  the  third 
harmonic  was  present  at  similar  rotor  speeds  (  «  1 1>~ 5  rpm)  with  a  stator  stagger 
angle  of  5". 2  degrees.  However,  its  level  in  that  case  was  not  strong  enough  to 
analyte  with  any  degree  of  accuracy. 

In  Ref.  1,  interference  noise  data  were  presented  for  the  same  rotor- 
stator  configuration  in  the  annular  cascade  as  used  in  the  current  test  program. 
The  maximum  rotor  speed  attained  in  those  past  tests  was  1450  rpm.  Only  the 
fourth  and  fifth  harmonics  of  blade  passage  frequency  propagated  at  these  low 
rotor  speeds.  However,  since  the  annular  cascade  had  been  disassembled,  modi¬ 
fied  to  increase  available  power,  and  reassembled  between  tests,  it  is  worth 
mentioning  that  the  current  test  results  agree  quite  well  with  the  past  test 
results  in  the  region  where  the  two  sets  of  data  overlap. 

It  is  evident  from  the  results  that  the  modified  annular  cascade 
attained  rotor  speeds  high  enough  to  excite  the  fundamental  as  well  as  higher 
harmonics  of  blade  passage  frequency.  Thus  the  experimental  acoustic  study 
succeeded  in  its  goal  of  providing  data  at  low  harmonics  of  blade  passage 
frequency  for  comparison  with  the  rotor-stator  interaction  noise  theory.  Com¬ 
parisons  between  this  theory  and  the  experimental  results  are  presented  in  the 
following  sub-section. 
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C.  COMPARISONS  WITH  THEORY 

I .  Theoretical  Model  And  Input  Data 

The  details  of  the  rotor-stator  interaction  noise  theory  were 
presented  in  Ref.  1,  along  with  the  earlier  comparisons  to  data  at  the  fourth 
and  fifth  harmonics  of  blade  passage  frequency.  The  acoustic  model  assumes 
that  the  rows  are  housed  in  an  infinitely  long,  hard-walled  annular  duct  of 
constant  hub/tip  ratio.  Each  row  is  represented  by  a  three-dimensional  array 
of  radially-oriented  line  dipoles,  whose  strength  is  fluctuating  in  time.  The 
resulting  acoustic  field  is  represented  as  a  superposition  of  propagating  duct 
modes  of  the  type  studied  by  Tyler  and  Sofrin.15  Reflections  from  neighboring 
blade  rows  or  terminations  of  the  duct  are  thus  neglected. 

The  strength  of  the  fluctuating  dipoles  is  determined  from  the 

rotor-stator  interaction  analysis  of  Osborne,  which  is  an  extension  to  com- 
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pressible  flow  of  the  classic  incompressible  model  of  Kemp  and  Sears. 

This  is  a  two-dimensional  model  which  assumes  that  the  unsteady  aerodynamic 
response  of  both  the  rotor  blades  and  stator  vanes  can  be  calculated  as  though 
they  operate  as  isolated  airfoils.36  The  results  of  Osborne's  analysis  are 
closed  form  expressions  for  the  fluctuating  lift  produced  by  each  of  three 
interaction  mechanisms:  unsteady  lift  on  the  rotor  produced  by  the  steady 
field  of  the  stator,  unsteady  lift  on  the  stator  produced  by  the  steady  field 
of  the  rotor,  and  unsteady  lift  on  the  stator  produced  by  the  impingement  of 
the  rotor  viscous  wakes.  The  first  two  are  usually  termed  potential  inter¬ 
actions,  and  the  last  a  viscous  interaction.  These  fluctuating  lift  predic¬ 
tions,  since  they  are  two-dimensional,  are  input  on  a  strip  theory  basis  to 
the  three-dimensional  acoustic  model  described  above. 

The  overall  theory  that  results  can  predict  the  total  acoustic 
energy  flux  away  from  the  stage  in  either  direction,  the  Sound  Pressure  Level 
(SPL)  at  any  point  in  the  duct,  and  the  individual  mode  amplitudes.  The  inputs 
required  are  the  steady  loading  on  both  rows,  which  is  used  in  the  potential 
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interaction  calculations ,  and  the  drag  coefficient  of  the  upstream  rotor  blade 
which  is  used  to  infer  the  wake  velocity  defect  for  the  viscous  interaction. 

As  part  of  a  larger  study  of  rotating  stall  phenomena, J  data  were  obtained 
on  the  circumferentially-averaged  turning  performance  and  total  pressure  loss 
across  the  rotor  operating  alone.  Viscous  losses  across  the  stator  were  as¬ 
sumed  negligible,  and  its  turning  performance  estimated  from  the  manufacturer' 
design  data.  The  earlier  comparisons  used  this  information  to  calculate  ap¬ 
proximate  values  for  the  inputs  required  by  the  theory,  using  Eqns.  (361 -(431 
of  Ref.  1. 


The  same  procedure  as  in  Ref.  1  was  used  in  the  latest  series 

of  comparisons,  except  that  the  reliability  of  our  input  dat3  was  improved 

in  several  respects.  The  main  improvement  is  that  we  now  have  data  on  the 

total  pressure  loss  across  the  rotor  in  the  presence  of  the  stator,38  whereas 

before  our  only  data  was  for  the  isolated  rotor.  As  a  result, -the  predicted 

Cp  of  the  rotor  blades  now  lies  in  the  range  0.05  -  0.0“,  which  is  much  more 

reasonable  than  the  value  0.16  which  had  been  predicted  from  the  isolated 

rotor  data.1  We  also  now  have  data  on  the  total  pressure  loss  across  the 

38 

stator  operating  as  part  of  the  st3ge,  whereas  before  no  data  on  the  stator 
were  available.  Although  not  needed  for  the  viscous  interaction  calculation, 
this  loss  does  play  an  indirect  role  in  estimating  the  stator  steady  lift  (see 
Eqn.  (42)  of  Ref.  1).  The  estimates  of  stator  turning  performance  are  now 
based  on  extrapolations  of  turning  data  measured  on  the  isolated  stator  row- 
in  the  annular  cascade,  whereas  in  the  earlier  study  manufacturers  design 
data  were  used  to  estimate  stator  turning.1  As  before,  we  had  to  rely  on 
isolated  rotor  turning  performance  due  to  the  difficulties  involved  in  trying 
to  instrument  the  small  gap  between  rotor  and  stator. 

2.  Discussion 


The  present  calculations  once  again  indicate  that  it  is  the 
rotor  viscous  wakes  impinging  on  the  stator  which  is  the  dominant  noise  genera 
tion  mechanism.  The  SPL  predictions  for  the  mean  stator  stagger  angle  of  28.2 
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are  compared  with  the  newly  acquired  data  for  the  first  through  fifth  harmonics 
of  blade  passage  frequency  in  Figs.  30a-d.  The  gaps  in  the  experimental  data 
are  those  regions  in  which  the  signal  was  close  enough  to  the  background 
level  that  no  reliable  level  could  be  inferred,  as  discussed  in  Section  III-B. 
In  each  figure,  the  rotor  rpm  at  which  the  various  duct  modes  first  propagate 
is  indicated  by  a  right-facing  arrow  with  a  pair  of  integers  which  denote  the 
azimuthal  and  radial  mode  numbers,  respectively.  It  should  be  recalled  that 
right  at  the  cut-off  rpm  the  mode  amplitude  is  theoretically  infinite,  since 
no  attempt  was  made  to  include  viscous  or  nonlinear  effects  in  the  model. ^ 

Thus  not  much  significance  should  be  attached  to  the  height  of  the  theoretical 
peaks  in  this  vicinity,  as  they  are  determined  primarily  by  how  closely  one 
wishes  to  approach  the  cut-off  rpm.  Generally,  the  calculations  were  made  at 
increments  of  50  rpm,  with  closer  spacing  near  the  cut-offs  to  better  exhibit 
where  the  modes  first  begin  to  propagate. 

Once  again,  the  highest  harmonics  are  underestimated  as  shown 
by  Figs.  30c,  d,  although  the  agreement  is  not  bad  at  the  higher  rpms  in  these 
two  cases.  Mote  that  this  is  because  the  theory  shows  a  definite  increase 
in  SPL  with  rotor  speed,  as  one  would  expect  intuitively,  whereas  the  average 
level  of  the  fourth  harmonic  data  changes  little,  and  the  fifth  harmonic  data 
show  only  a  slight  rise.  In  general,  one  can  also  reasonably  correlate  the 
location  of  the  theoretical  peaks  near  cut-off  with  those  in  the  experiments. 
However,  the  data  do  show  much  more  of  a  peak/valley  structure  than  does  the 
theory.  Mo  firm  explanation  is  apparent  for  this,  although  it  may  be  the 
result  of  the  reflection/transmission  of  the  acoustic  signals  by  the  upstream 
rotor,  which  is  not  accounted  for  in  the  theory. 

Curiously,  the  comparisons  at  the  lower  harmonics,  Figs.  30a, 
b,  show  a  reversal  in  the  levels  of  theory  and  experiment.  The  predictions 
are  now  too  high  over  the  entire  range  in  which  data  could  be  taken,  which 
was  limited  again  by  power  constraints.  The  cut-on  of  the  second  harmonic 
is  reasonably  well  predicted;  little  can  be  said  of  the  third  harmonic  because 
of  the  very  limited  range  in  which  it  could  be  distinguished  in  the  third- 
octave  band  spectrum.  A  notable  feature  of  the  comparisons  appears  at  the 
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that  the  level  of  the  interaction  noise  at  these  harmonics,  if  excited,  will  he 
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fundamental  blade  passage  frequency;  not  only  are  the  predictions  too  high, 
but  the  data  show  that  the  fundamental  is  propagating  in  the  duct  well  before 
the  theory  would  predict,  at  least  on  the  basis  of  rotor-stator  interaction. 

From  the  present  data  one  cannot  tell  whether  the  fundamental  first  cuts  on 
at  about  1490  rpm,  or  whether  it  had  been  present  before,  but  below  the  back¬ 
ground  level  (see  Fig.  26).  In  any  event,  this  strongly  suggests  that  another 
noise  generation  mechanism  besides  rotor-stator  interaction  may  be  influencing 
the  data.  This  will  be  discussed  further  below. 

The  comparisons  of  theory  with  experimental  data  at  the  mean 
stator  stagger  angle  of  37.2°  are  shown  in  Figs.  31a-d.  Qualitatively,  the 
comparisons  exhibit  the  same  features  as  those  for  the  other  stator  setting. 

The  predictions  for  the  fourth  and  fifth  harmonics  are  too  low,  though  the 
agreement  at  the  higher  rpms  is  not  bad  in  the  case  of  the  fourth.  The  lower 
harmonics  are  overestimated,  though  to  a  lesser  degree  than  before.  Again, 
the  fundamental  is  detectable  well  before  the  theory  indicates  it  should 
propagate. 

As  noted  in  connection  with  the  comparisons  done  under  the 
previous  program,^  one  plausible  explanation  for  the  poor  agreement  at  the 
higher  harmonics  is  the  assumed  shape  for  the  viscous  wake  profiles.  At 
present,  the  model  uses  a  profile  correlation  which  is  based  on  the  isolated 
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airtoil  data  ot  Silverstein,  et  al.  '  It  was  demonstrated  in  Rot.  1  that 
alternate  correlations  of  the  same  data  would  produce  significant  changes  in 
the  theoretical  SPL  levels  of  the  higher  harmonics.  This  is  basically  because 
the  high  harmonic  content  of  the  wake  profiles,  and  hence  the  noise,  is  deter¬ 
mined  to  a  great  extent  by  the  rate  at  which  the  velocity  defect  approaches 
cero  near  the  edges  of  the  wake,  and  this  is  difficult  to  define  experimentally. 
Thus,  it  may  be  that  reliable  predictions  of  viscous  interaction  noise  in 
the  higher  harmonics  will  have  to  await  the  development  of  means  to  predict 
the  wake  structure  with  more  precision. 

In  a  practical  sense,  the  inability  to  predict  higher  harmonic 
levels  is  of  secondary  interest.  It  is  the  fundamental,  and  perhaps  second 


harmonic,  which  is  most  directly  responsible  for  the  noise  problem  at  the  higher 
speeds  at  which  aircraft  engines  operate.  Here  the  theory  consistently  over¬ 
estimates  the  noise  levels.  This  cannot  be  explained  on  the  basis  of  inade¬ 
quate  knowledge  of  the  wake  profile  structure;  as  demonstrated  in  Ref.  1,  the 
lowest  harmonics  are  insensitive  to  the  particular  choice  of  profile  correla¬ 
tion.  They  are,  however,  sensitive  to  the  prediction  of  the  wake  centerline 
velocity  defect,  i.e.,  the  height  of  the  wake  profile.  This  also  was  correlated 
on  the  basis  of  the  isolated  airfoil  data  in  Ref.  39.  More  recent  experimental 
investigations  have  been  accumulating  data  on  the  wake  structure  of  both  two- 
dimensional  cascades  and  three-dimensional  rotors.  The  cascade  data 

of  Ref.  40  indicate  that  the  decay  of  the  wake  velocity  defect  with  distance 
downstream  of  the  trailing  edge  is  slower  than  would  be  found  behind  an  iso¬ 
lated  airfoil.  On  the  other  hand,  Kerrebrock  et  alJ1  observed  a  faster  rate 
of  decay  in  the  cascade  data  when  compared  with  the  isolated  airfoil  case. 
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The  studies  done  on  three-dimensional  rotor  wakes  *’  thus  far  have  shown 
that  they  decay  more  rapidly  than  the  wakes  of  either  cascades  or  isolated 
airfoils. 

Clearly,  there  is  ample  evidence  to  suggest  that  the  present 
theory  may  be  inaccurately  modeling  the  rotor  wake  structure.  Improved  cor¬ 
relations  with  the  present  acoustic  data  would  result  if  the  centerline  ve¬ 
locity  defect  was  weaker  than  predicted,  and  the  wake  profile  shape  was  sharper 
and  narrower  than  that  assumed.  However,  it  is  too  early  to  say  how  generally 
applicable  the  data  of  Refs.  40-43  are  to  other  configurations.  In  particular, 
the  effects  of  such  parameters  as  solidity,  steady  loading,  and  blade  profile 
all  need  to  be  quantified.  Lacking  this,  better  knowledge  of  the  wake  struc¬ 
ture  in  the  present  annular  cascade  facility  would  require  an  extensive  series 
of  measurements  which  were  felt  to  be  outside  the  scope  of  the  present  pro¬ 
gram. 

Another  possible  source  of  error  in  the  theory  is  the  unsteady 
aerodynamic  model.  As  pointed  out  earlier,  Osborne's*''’  model  does  include 
compressibility  effects  provided  the  Mach  number  and/or  reduced  frequency  are 
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were  measured  on  the  outer  casing  upstream  of  the  rotor.  In  taking  these  data, 
the  fan  system  downstream  of  the  annular  cascade  was  turned  off  and  the  fan  was 

SO 


* 


not  too  high.  However,  when  applied  at  the  low  speeds  of  the  present  experi- 
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ments,  it  reduces  to  the  incompressible  theory  ot  Kemp  and  Sears.'  ’  The 
flow  is  assumed  to  be  two-dimensional,  with  each  blade  or  vane  behaving  as 
an  isolated  airfoil.  For  the  high  hub/tip  ratio  of  the  present  experiments, 
the  two-dimensional  assumption  is  not  unreasonable.  However,  the  solidity 
ratio,  defined  as  the  ratio  of  chord  length  to  tangential  blade  spacing,  is 
approximately  0.S5  for  the  stator  at  mid-annulus.  Thus  the  effects  of  inter¬ 
ference  from  neighboring  stator  vanes  may  not  be  negligible.  Both  of  these 
questions,  regarding  the  importance  of  three-dimensionality  and  blade-to- 
blade  interference,  could  be  answered  by  the  type  of  analysis  reported  in 
Section  II,  which  is  still  under  development.  The  best  means  of  assessing 
the  aerodynamic  model  would  be  to  make  unsteady  load  measurements  on  the 
stator  vanes  and  see  how  they  compared  with  predictions.  This  was  also  felt 
to  be  outside  the  scope  of  the  present  work. 

Finally,  there  is  the  question  of  why  the  fundamental  blade 
passage  frequency  cuts  on  earlier  in  the  experiments  than  the  theory  would 
indicate  (Figs.  30a  and  31a).  This  strongly  suggests  that  noise  is  being 
generated  through  mechanisms  other  than  rotor-stator  interaction.  The  most 
likely  possibility  is  that  the  rotor  itself  is  interacting  with  a  distortion 
in  the  inlet  flow,  perhaps  the  wakes  produced  by  the  upstream  support  struts 
shown  in  Fig.  14.  Since  the  struts  are  streamlined  and  located  about  a  foot 
and  a  half  ahead  of  the  rotor,  one  would  expect  that  their  wakes  would  be 
sufficiently  diffuse  by  the  time  they  reach  it  so  as  to  have  neglibible  in¬ 
fluence.  Azimuthal  surveys  of  the  flow  ahead  of  the  rotor  have  been  performed, 
but  onlv  in  the  presence  of  distortion  screens  which  were  used  to  study  the 
influence  of  distortion  on  rotating  stall.'1  Any  strut  generated  distortion 
would  have  been  masked  by  the  presence  of  the  screens. 

Thus  it  seemed  the  most  expeditious  way  to  test  this  hypothe¬ 
sis  was  to  remove  the  stator  and  see  whether  the  isolated  rotor  still  produced 
significant  tone  noise.  (In  the  absence  of  inflow  distortion,  it  should  not 
radiate  any  noise  at  subsonic  tip  speeds.)  When  the  rig  was  disassembled  for 


mechanism  to  fire  on  the  audio  oscillator  signal  and  then  adjusting  the  oscillator 
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this  purpose,  it  was  discovered  that  the  chain  drive  used  to  drive  the  rotor 
had  deteriorated  and  loosened  considerably.  Accordingly  it  was  replaced,  but 
the  new  drive  generated  sufficiently  higher  background  noise  levels  than  the 
old  that  further  experiments  on  isolated  rotor  noise  were  not  possible. 

D.  CONCLUSIONS 

Our  earlier  comparisons  of  the  rotor-stator  interaction  noise  theory 
with  experimental  data  taken  in  Calspan’s  annular  cascade  facility  had  shown 
that  the  theory  consistently  underestimated  STL  at  the  fourth  and  fifth 
harmonics  of  blade  passage  frequency.''  This  was  felt  to  be  a  rather  severe 
test  of  the  theory,  as  the  higher  harmonics  were  shown  to  be  quite  sensitive 
to  details  of  the  rotor  blade  wake  profiles.  However,  power  limitations  in 
the  rig  at  that  time  prevented  operating  at  high  enough  rpm's  to  allow  the 
lower  harmonics  to  propagate.  The  main  goal  of  this  portion  of  the  current 
contract  was  to  provide,  with  a  newly  re-powered  rig,  comparisons  between 
theory  and  data  at  the  fundamental  and  first  few  harmonics.  These  comparisons 
have  been  made,  and  show  that  the  theory  consistently  overestimates  the  noise 
in  the  fundamental  and  second  harmonic;  experimental  data  at  the  third  har¬ 
monic  was  also  lower  than  the  theory,  but  was  limited  to  too  few  points  for 
any  general  conclusion  to  be  drawn.  Interaction  of  the  stator  with  the  viscous 
wakes  shed  by  the  rotor  is  predicted  by  the  model  to  be  the  predominant  noise 
generation  mechanism;  potential  interactions  between  blade  rows  were  found 
to  be  negligible. 

The  theoretical  model  was  required  to  predict  the  viscous  interaction 
noise  based  only  on  turning  and  loss  performance  data,  since  this  is  the  type 
of  information  the  designer  is  likely  to  have  available.  The  overall  model 
is  thus  a  composite  of  three  major  elements:  a  model  of  the  viscous  wake 
structure;  an  aerodynamic  model  to  predict  the  unsteady  response  of  the 
downstream  row  to  the  wakes;  an  acoustic  model  to  predict  the  resulting  radi¬ 
ation.  At  present  we  know  with  certainty  only  that  the  end  results  do  not 
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ithin  the  boundary  layer  on  the  outer  wall  and  by  the  free-wheeling  fan  down- 
tream  of  the  annular  cascade  test  section.  Figs,  l'b  and  l'c  show  the  sharp 

S2 


correlate  well  with  the  acoustic  measurements,  h’e  strongly  suspect  that  the 
weakest  links  in  the  model  are  the  viscous  wake  modeling,  and  perhaps  the 
unsteady  aerodynamic  model,  for  the  reasons  discussed  above.  Conclusive 
answers  to  these  questions  can  be  obtained  only  by  performing  extensive 
measurements  of  both  the  rotor  wake  structure  and  the  unsteady  loads  produced 
on  the  stator  vanes;  however,  such  a  program  was  well  outside  the  scope  of 
present  effort. 

Another  surprising  feature  in  the  comparisons  is  that  the  fundamenta 
blade  passage  frequency  first  appears  at  a  lower  rpm  than  predicted  by  the 
theory.  This  suggests  that  a  noise  generation  mechanism  other  than  rotor- 
stator  interaction  may  have  heavily  influenced  the  data.  One  possibility  is 
that  the  rotor  itself  may  be  interacting  with  a  distortion  in  the  flow,  such 
as  wakes  from  the  upstream  support  struts,  or  perhaps  a  distortion  in  the 
inlet  flow  itself.  N'o  surveys  have  been  made  as  yet  of  the  flow  uniformity 
ahead  of  the  rotor,  and  attempts  at  measuring  the  noise  output  of  the  rotor 
with  the  stator  removed  were  thwarted  by  high  background  levels  from  a  newly 
installed  drive  mechanism.  It  is  felt  that  one  or  both  of  these  measurements 
should  be  made  before  conducting  any  further  experiments  on  rotor-stator 
interaction  in  the  annular  cascade. 


- tltMuvuv* .  .*\u  example  or  tne  procedure 

xs  shown  in  Fig.  13a.  In  this  case,  the  relative  background 
level  for  the  third-octave  band  centered  at  1250  H;  is  estimated 


SECTION  IV 
SUMMARY 

In  Section  I,  a  method  is  presented  for  computing  the  steady  loading 
on  a  three-dimensional  annular  blade  row  in  compressible  flow,  given  the  blade 
camber  line.  A  kernel-function  procedure  is  used  to  solve  the  linearised  in¬ 
tegral  equation  relating  the  unknown  loading  to  the  camber  line  slope.  The 
greatest  advantage  of  the  solution  procedure  presented  here  is  that  it  allows 
analytical  expressions  to  be  derived  for  both  the  spanwise  and  chordwise  in¬ 
tegrals  which  are  required.  It  also  allows  the  various  aerodynamic  and  per¬ 
formance  parameters  of  the  rotor  to  be  expressed  as  simple  algebraic  functions 
of  the  loading  expansion  coefficients. 

Comparisons  are  presented  between  the  predictions  of  the  present 
tneory  and  those  of  two-dimensional  strip  theory,  as  well  as  with  the  results 
of  an  inverse  three-dimensional  analysis.  The  three-dimensional  and  strip 
theories  show  consistently  good  agreement  at  mid-annulus,  but  the  three- 
dimensional  theory  predicts  loadings  which  are  significantly  higher  near  the 
hub  and  lower  near  the  tip  than  the  strip  theory,  thus  giving  much  lower 
radial  variations.  Such  behavior  can  be  attributed  to  the  induced  angle  of 
attack  distribution  generated  by  the  trailing  vortex  wakes,  which  is  not 
accounted  for  in  strip  theory;  this  smoothing  appears  to  be  the  principle 
effect  of  three-dimensionality,  at  least  for  subsonic  relative  Mach  numbers. 

When  coupled  with  the  analogous  treatment  of  the  disturbance  flow- 
field  produced  by  blade  thickness,  the  loading  analysis  presented  here  should 
provide  an  efficient  means  of  studying  three-dimensional  effects  in  turbo¬ 
machinery  rows.  The  program  should  prove  particularly  advantageous  for 
evaluating  alternate  blade  camber  lines  at  a  given  operating  condition;  thus, 
it  could  serve  to  narrow  down  the  number  of  candidate  blade  profiles  before 
going  to  a  finite-difference  solution  of  the  full  nonlinear  equations.  It 
could  also  prove  useful  in  providing  better  initial  conditions  to  a  difference 
scheme,  allowing  the  latter  to  converge  in  fewer  iterations. 
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In  Section  II,  the  linearized  analysis  is  extended  to  the  study  of 
unsteady  flow  through  a  rotor,  such  as  might  be  produced  by  inflow  distortion 
or  blade  vibrations.  An  integral  equation  is  derived  which  relates  the  un¬ 
known  loading  to  the  prescribed  variations  in  normal  velocity  at  the  blade 
surfaces.  Terms  can  be  identified  in  the  integral  equation  which  represent 
the  blade  wakes.  These  include  elements  of  both  trailing  vorticity,  produced 
by  spanwise  loading  gradients,  and  shed  vorticity,  which  results  from  the 
temporal  fluctuations  in  loading.  Also  included  are  terms  which  represent 
waves  propagating  unattenuated  away  from  the  blade  row  in  both  directions. 

This  acoustic  radiation  field  would  be  predicted  as  a  by-product  of  the  solu¬ 
tion,  with  regard  to  both  total  acoustic  energy  flux,  and  its  distribution 
among  the  various  propagating  duct  modes.  Though  somewhat  more  complicated 
than  the  analogous  steady  flow  equation,  the  form  of  the  required  integrations 
suggests  that  it  can  be  solved  by  a  suitable  extension  of  the  procedure  de¬ 
veloped  in  Section  I. 

Section  III  presents  comparisons  between  a  theoretical  model  of 
rotor-stator  interaction  noise,  developed  under  a  previous  program,  and  a 
current  series  of  acoustic  measurements  made  in  an  annular  cascade  facility. 
Inputs  required  by  the  theory  are  the  steady  loading  on  both  rows,  which  is 
used  to  predict  the  potential  interactions ,  and  the  drag  coefficient  of  the 
upstream  row,  which  is  used  to  infer  the  wake  velocity  defect  for  the  viscous 
interaction.  These  were  inferred  from  experimental  measurements  of  the  cir¬ 
cumferentially-averaged  turning  performance  and  total  pressure  losses, 
respectively ,  across  the  blade  rows. 

Theoretical  and  experimental  Sound  Pressure  Levels  at  one  point  in 
the  duct  upstream  of  a  rotor-stator  pair  were  compared  from  the  fundamental 
blade  passage  frequency  up  through  the  fifth  harmonic.  The  calculations  in¬ 
dicate  that  the  viscous  wake  interaction  is  the  dominant  noise  generation 
mechanism.  The  theory  overestimates  the  noise  in  the  fundamental  and  second 
harmonic,  while  its  predictions  are  too  low  for  the  fourth  and  fifth.  The 
most  likely  sources  of  error  in  the  theory  are  felt  to  be  in  the  modeling  of 
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the  rotor  viscous  wakes  and  the  unsteady  aerodynamic  response  of  the  down¬ 
stream  stator  vanes;  the  wake  model  is  based  on  isolated  airfoil  data,  while 
the  aerodynamic  analysis  neglects  the  influence  of  neighboring  vanes  in  the 
stator.  Detailed  wake  survey  and  unsteady  load  measurements  against  which 
to  test  this  hypothesis  were  outside  the  scope  of  the  present  effort.  It 
also  appeared  that  the  experimental  data  may  have  been  heavily  influenced 
by  noise  generation  mechanisms  other  than  those  considered  in  the  rotor- 
stator  interaction  theory,  including  the  possibility  of  the  rotor  interacting 
with  the  wakes  shed  by  upstream  support  struts. 


Appendix  A. 


EQUIVALENCE  OF  THE  PRESSURE  DIPOLE  .AND  HORSESHOE  VORTEX 
REPRESENTATIONS  OF  THE  STEADY  LOADING  ON  AN  ANNULAR  BLADE  ROW. 


The  pressure  dipole  representation  of  the  rotor  used  in  the  steady 


lifting-surface  calculations  can  be  shown  to  be  equivalent  to  the  vortex 
_ _ _ _ _ i _  Ai. _  _  i  » .  ^ _  S  —  7  >  i  .  * 


representation  used  by  Okurounmu  and  McCune  .  In  order  to  demonstrate  this, 
the  present  dipole  representation  is  reinterpreted  as  a  vortex  representation, 
to  which  it  must  be  equivalent  if  both  are  valid.  The  resulting  vortex  repre¬ 
sentation  is  then  shown  to  be  the  same  as  that  of  Okurounmu  and  McCune. 


First,  the  perturbation  velocity  potential  of  a  lifting  rotor  is 
written  in  terms  of  the  perturbation  pressure,  -p  ,  as  follows: 
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Writing  the  perturbation  pressure  as  the  superposition  of  dipole  solutions 
and  then  interchanging  the  order  of  integration  gives 
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where  Sg  means  all  blade  surfaces  . 
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Introducing  the  bound  vorticity,  #  ,  which  is  related  to  the  blade  loading 

by  T{  *  -  &  /pM  I),,  , 
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The  sign  of  the  vorticity  is  defined  opposite  to  that  in  the  main  body  of  the 
report  to  be  consistent  with  Okurounmu  and  McCune.  Now  it  can  be  shown  that 


J  f>0  (Z\ Z,  ,  r  ;  ,  r#,  2J.  )  d  £  = 


~P  (  %  >  ^  t  ^  4  -  >  ^  a)  ^  a 


Also,  the  perturbation  pressure  and  potential  fields  both  must  satisfy  the 
same  partial  differential  equation  and  boundary  conditions  (see  Eqns .  (55) 
and  (56)  of  Ref.  1).  Thus,  it  follows  that  the  solution  for  a  pressure  dipole, 
p0  ,  is  mathematically  the  same  as  that  for  a  potential  doublet,  , 
although  the  two  represent  physically  different  flows.  Analogous  to  the 
situation  in  wing  theory  (Ref.  16,  Chapt.  5),  we  can  then  derive  the  potential 
due  to  an  elementary  horsehoe  vortex  as: 
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Then  Eqn.  (A-4)  becomes 
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This  equation  expresses  the  potential  due  to  a  lifting  rotor  as  the  super¬ 
position  of  elementary  horseshoe  vortices.  Each  vortex  has  an  infinitesimally 
long,  radially  oriented  bound  element  located  at  z0  ,  r„  and  a  pair  of 
trailing  vortex  filaments  lying  along  the  helical  undisturbed  streamlines. 


In  Ref.  5,  Okurounmu  and  McCune  developed  the  solution  for  equally 
spaced,  radially  oriented  vortex  lines  (and  their  associated  trailing  vortex 
wakes)  which  span  the  annulus.  To  show  that  their  solution  is  equivalent  to 
that  based  on  a  dipole  formulation,  we  proceed  as  follows.  First,  we  sub¬ 
stitute  our  pressure  dipole  solution,  Eqn .  (91)  of  Ref.  1,  for  0  in 
Eqn.  (A-5) .  The  required  integration  has  already  been  performed  in  Ref.  1, 
and  is  quoted  there  as  Eqn.  (100a).  This  elementary  horseshoe  vortex 
solution  can  then  be  used  in  (A-6)  to  construct  the  potential  field  due  to 
B  radially-oriented  line  vortices  whose  strength,  f  ,  varies  with  span: 


Carrying  out  the  required  substitutions  and  the  sum  over  ^  ,  we  obtain 


This  solution  will  now  be  shown  to  be  equivalent  to  that  derived  by 
Okurounmu  and  McCune^  for  the  same  problem,  though  they  worked  from  a  different 
point  of  view.  When  notational  differences  are  accounted  for  the  nn  =•  o  terms 
in  the  above  expression  are  easily  shown  to  be  identical  to  those  in  Eqn.  (8) 
of  Ref.  5.  However,  the  key  step  in  demonstrating  that  the  two  solutions  are 
equivalent  is  in  showing  that  the  terms  representing  the  trailing  vortex  wakes 
are  the  same.  These  are  the  mto  terms  in  £Pp  which  appear  in  the  second 
and  third  lines  of  Eqn.  (A-S) ;  they  will  be  denoted  by  for  convenience. 

Using  Eqn.  (2a),  the  wake  potential  may  be  written: 
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We  now  employ  the  addition  and  subtraction  manipulation  used  in  Section  IV-C 
of  Ref.  1  to  demonstrate  that  the  expressions  for  the  velocity  components 
contained  the  correct  discontinuities.  Then  becomes 
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The  sum  over  k  in  the  first  line  is  the  Fourier-Bessel 
1  <~r'  ;  the  sum  over  rv\  yields  the  sawtooth  function  X-  defi 


sset  expansion  tor 
defined  by: 
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Note  that  <->  ^  has  a  discontinuity  of  (  —  -g—  )  as  one  crosses  a  blade 


surface,  ^  >  in  the  positive  direction.  By  substituting  from 

Bessel's  equation,  the  second  line  of  (A-10)  can  be  written  in  terms  of  de¬ 
rivatives  of  the  series  S„  (r,  r#  )  defined  by: 
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The  wake  potential  then  reads 
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The  series  defined  in  (A- 1 2)  has  been  summed  by  Salaun1->  in  terms 
of  modified  Bessel  functions  using  the  following  technique.  If  we  differen¬ 
tiate  the  series  twice  with  respect  to  <P  ,  and  use  the  Fourier-Bessel  expan¬ 
sion  for  S  ( <T-<r  ) , 
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then  the  following  differential  equation  is  obtained  for  : 
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LA- 15) 


where  the  notation  d"^  =  noor/l/  has  been  introduced.  (Note  that  the  addition 
and  subtraction  manipulation  used  in  Section  IV-C  of  Ref.  1  is  equivalent  to 
substituting  for  from  Eqn .  (A- 15)).  This  equation  is  solved  subject  to 
the  hard-wall  boundary  conditions,  i.e.  that  d.  /  oi<TV,  vanish  at  the 
inner  and  outer  duct  walls  <X^,  =  r,u)lZ  '  ^  aJtd  CT,-,t  =  n  ujr~T  ’  ^  •  Moreover, 

is  required  to  be  continuous  at  <3%,  *  crn<)  and  cLS^/dv^  must  satisfy  the 
following  jump  condition  there: 
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This  was  derived  by  integrating  Eqn.  (A- 15)  from  £  to  J~no  +£ 

The  modified  Bessel  functions  InC<rn)  and  Xn(<rn)  are  linearly 
independent  solutions  of  the  homogeneous  equation  for  S  „  •  For  <Tn  <  crn 
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The  coefficients  C(  thru  C ^  are  chosen  so  as  to  satisfy  the  boundary 
conditions  cited  above,  and  the  results  are: 
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and  primes  denote  differentiation  with  respect  to  the  argument. 

Returning  to  Eqn .  (A-13),  if  we  integrate  by  parts  once  in  the 
second  term,  make  use  of  the  fact  that  dS^&/dra  vanishes  at  the  hub  and 
tip,  and  then  substitute  for  5  me  trom  the  results  given  in  (A- 17)  and 
(A- IS),  we  obtain  for  the  wake  potential: 
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After  reorganizing  the  ranges  of  the  integrals  above,  the  expression  in  braces 

can  be  shown  to  be  >  where  is  the  wake  function  defined 

(3  i +on 

in  Eqn.  (5)  of  Ref.  5.  (Okurounmu  and  McCune  placed  their  blades  at  ^  - 

2  TT  -l  m 

rather  than  — ,  which  results  in  an  extra  factor  of  (-1)  in  their 
o 

result.)  The  above  representation  for  the  wake  potential  is  thus  equivalent 
to  that  given  by  Okurounmu  and  McCune5  in  Eqn.  (4). 


It  remains  to  be  shown  that  the  last  four  lines  in  Eqn.  (A-S) , 
which  contain  the  exponential  factors  in  (  2  -  ) ,  are  the  same  as  the 

corresponding  ones  in  the  line  vortex  solution  of  Okurounmu  and  McCune.  The 
equivalence  of  the  m  -o  terms  is  straightforward,  as  noted  earlier.  For 

,  the  factors  in  the  braces  in  Eqn.  (A-8)  are  first  rearranged,  using 
the  same  addition  and  subtraction  manipulation  referred  to  above.  Then,  the 
following  two  radial  integrals  which  arise, 

rc<ra)  R^3k  (<r4 )  d  <ra 


_  [  2  z  (  Krv,0k  \z' 


f  /  Km»*  Ya  . 

L(  no 6  J 

{  u  ) 

d<r» 


can  be  identified  as  the  coefficients  of  (7~)  in  the  Fourier-Bessel 

expansions  of  the  bound  vorticity,  P(cT)  ,  and  the  wake  functions,  X^S> 
respectively.  The  latter  follows  from  the  above  demonstration  that  the  A 
series  in  the  second  ten.,  of  Eqn.  (A-10)  is  indeed  proportional  to  \yV,g 
Okurounmu  and  McCune  refer  to  the  coefficients  defined  by  the  above  integrals 
as  r’miiu  and  >  respectively.  Once  these  identifications  are  made, 

the  equivalence  between  these  terms  and  their  counterparts  in  Eqn.  (8)  of 
Ref.  5  is  evident. 

Thus,  we  have  been  able  to  show  that  the  singularity  solutions  used 
to  represent  a  lifting  rotor  in  the  present  work  can  be  put  in  the  same  form 
as  those  used  by  Okurounmu  and  McCune5.  On  this  basis,  we  have  concluded  that 
the  lifting-surface  formulations  based  on  the  pressure  dipole  and  vortex  rep¬ 
resentations,  though  in  markedly  different  forms,  are  indeed  equivalent. 
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Appendix  B 

EVALUATION  OF  RADIAL  INTEGRALS 


This  appendix  presents  the  analytical  evaluation  of  the  set  of  radi¬ 
al  integrals  defined  by: 


j  r ' 

n  *  *  /  ^  \  k 


where  ,  k  and  n»rm8  are  all  integers.  The  are  linear  combinations  of 
the  J„  and  yn  Bessel  functions: 


^ )  ~  d  2s) 


where  c5‘n)<  and  Mnk  are  phase  and  normalitation  factors,  respectively,  which 
vary  from  mode  to  mode,  but  are  independent  of  <T  .  These  functions  are 
discussed  at  length  in  Refs.  2-4,  and  15.  The  procedure  used  in  evaluating 
the  eigenvalues  and  eigenfunctions  for  large  values  of  n  is  described  in 


The  above  integral  can  be  evaluated  by  using  the  following  indefinite 


integral  given  by  Watson:' 


j  f  )  d  f  =  ( -J  ■> ■fl-l)  f  )  S^_,  ^  (D  “  <  -n~<  (?)  )  (B-5) 

where  2n  represents  either  J0  or  Yn  and  the  are  Lommel's  functions. 

Substitution  of  Eqns.  (B-2)  and  (B-3)  into  (B-l)  yields: 

I*k  *  — -  -  (j -')  2  n«(  t*,*  -l(K~K<r') 

f 

[J«-,  (  ■k’rtK  <T  )  “  Sni<Yn.1  c  fu 

-  - -z -  f 

1  J  Ja. 


This  can  be  simplified  by  making  use  of  the  Bessel  function  recursion  rela¬ 
tion, 

-  Op  -  fzjv 

where  a  prime  indicates  differentiation  with  respect  to  the  argument.  Also, 
by  definition: 


—  Z„k(*„k<r)  -  o  <r  •  A,  i 

o.<r 


This  gives 


S  >w./^kr)-  ~  S  {K^r)  j* (B-4) 

A nk  u  nk 3  J  _ A 


Watson"  also  gives  the  following  recursion  formula  satisfied  by  Lommel's 


function : 


aUd — f?)  -  f 


which  allows  Eqn.  (B-4)  to  be  put  in  its  simplest  form: 


f  <rZ«k(Knkr)  S' tnCK„k  <r) 


(B-5) 


Since  Rnl<  at  the  hub  and  tip  will  have  already  been  evaluated  as 
part  of  determining  Knk  ,  the  only  new  quantities  introduced  above  are  the 
Lommel  functions.  They  may  be  easily  evaluated  using  the  formulae  cited  by 

23 

Watson;"  since  these  are  all  given  for  vS^>ri  rather  than  its  derivative,  in 
practice  Eqn.  (B-4)  was  used  rather  than  (B-5).  The  calculated  values  were 
checked  by  recognizing  that  1^  as  defined  in  Eqn.  (B-l)  is  just  the  k  th  co¬ 
efficient  in  the  Fourier- Bessel  expansion  of  ’  ,  i.e. 


%nk(Knk  r') 


(B-6) 


Satisfaction  of  Eqn.  (B-6)  was  verified  for 


several  values  of  ^  ,  n  and  <p  . 
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Appendix  C 

EVALUATION  OF  AXIAL  INTEGRALS 


This  appendix  presents  the  analytical  evaluation  of  the  axial  in¬ 
tegrals  defined  by  Eqn.  (14).  C\f>  is  easily  evaluated  after  making  the 
substitution 


=  £<W  o. 

The  results  are: 


X  a  c*tn 


(C-l) 


rt  -  <p  +  a 

■  t  '  / 

2  to  /  4  • 

*  *  c 

. 

t  4  /2  4  ~  •**'"  C  i  - 2 )  <J>  /2(i  -2)  i  >2 


The  derivation  of  the  results  for  C.'m&k  and  5^ &k  is  rather  lengthy, 
although  the  procedure  is  straightforward  once  a  few  key  relations  are  given. 
Only  an  outline  of  the  principal  steps  is  given  below,  followed  by  the  final 
results : 

1)  The  cosines  and  sines  appearing  in  the  integrands  of  Eqn.  (14) 
are  written  as  either  the  real  or  imaginary  parts  of  complex  exponentials, 
and  the  region  of  integration  is  split  in  two,  according  to  whether  %  is 
greater  or  less  than  %  . 

2)  Use  is  made  again  of  the  substitution  in  Eqn.  (C-l);  this  reduces 
the  integrands  to  trigonometric  forms  involving  only  sines  and  cosines  of 
multiples  of  .  For  t  >  /  ,  set 

c  l  <Pt  M  2*  [  2*04  C  I  -  2  )  ipa  ~  2*04  L  J 

l  <p 

3)  Each  integrand  will  contain  an  exponential  of  the  form  £ 
where  H  is  in  general  complex.  It  takes  on  the  values  2  *  ia-cb  where 


SI 


The  exponential  is  expanded  in  a  simple  cosine  series  using  Eqn.  (9.6.34)  of 


Here  -£,,(2)  is  the  modified  Bessel  function  of  the  first  kind  of  integer  order 

n  .  N'oting  that  lln(.£)l  =  /  Jn  (Z )  I  ,  where  Jn(2>  is  the  regular  Bessel 

function,  and  that  decays  exponentially  with  n  once  o  ~  1 1[  (see, 

25 

e.g.,  Michels,  )  the  series  in  Eqn.  (C-4)  can  be  expected  to  converge  rapidly 
once  n  exceeds  this  value. 


4)  The  integrand  will  now  contain  a  series  of  terms,  each  of  which 
involves  at  worst  the  product  of  two  cosines.  Since  we  know  the  series  con¬ 
verge  as  discussed  in  3),  the  integration  over  <p0  may  be  carried  out  term-bv 
term. 


5)  The  series  of  integrated  terms  which  result  from  4)  will  have  dif¬ 
ferent  arguments  for  ,  depending  on  whether  the  series  arose  from  the 
integral  for  Xa  greater  or  less  than  %  .  However,  the  two  are  related  through 


where  the  asterisk  denotes  the  complex  conjugate.  This  allows  most  of  the 
terms  in  the  series  arising  from  the  two  integrals  to  be  combined  in  a  single 
series,  considerably  simplifying  the  results. 


The  final  results  are  given  below,  where  the  symbols  A c  and  denote 
eal  and  imaginary  parts,  respectively: 


i ^  3k  '  *"  ^ 


c^noX  Re.!,)  -  d  X  ( Za  -  1 ,  ) 


f  Z  \  -On-"’  <?  -  £  )  ( iJS-4  aX  O-W  0  Z  Jle  _  r  ^ 


♦  ^  c-o 

r>  •  f 


2  0*>  - 1 ) 


f  J ax  1  . 

<  .  >  C-«  0  X  tfCg.  I  _  +■ 

1  i^~~Jlaxi 


XeZa  t  ^ 

2~/L  b  Z  >*  Z. 

.3+**  n  O 

Lo^v.  (nri)  <£  "J 

n 

2  (n  n  O  j 

[_  OoA-  a  X 

"j  1 

>  i««io  X  Jy^,  A„ 

J  J 

i  r  \  . 

-  cjX  C  -  &  I,  j 

n  OOO 
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where  now  1  n  a  '■'<J3  °  *  °x  ^  -«• 

In  the  above  expressions  <5^  is  the  Kronecker  delta  symbol,  defined  as  one 

when  l  -  j  ,  and  zero  otherwise.  It  must  also  be  remembered  in  evaluating  the 

sums  over  o  that  JU~*  '^ISl  ^  =  a  .  The  special  case  m  =  o  can  be  obtained 
n  —  o  ^ 

from  the  above  results  by  setting  b-o  and  recognizing  that  this  implies 
the  Bessel  functions  will  be  purely  real. 

At  first  glance,  the  evaluation  of  Eqn.  (C-6)  might  appear  even 
more  time-consuming  than  a  more  straightforward  numerical  integration.  But, 
in  addition  to  the  obvious  advantage  that  they  are  "exact"  to  any  desired 
accuracy,  depending  on  how  many  terms  one  retains  in  the  sum  over  n  ,  there 
is  one  overriding  point  in  their  favor.  That  is,  whereas  a  numerical  integra¬ 
tion  scheme  must  repeat  every  step  of  the  calculation  for  all  combinations 
of  inn  ,  k  ,  i  and  2  ,  the  above  expressions  conveniently  separate  the  influence 
of  each  of  these  indices.  For  example,  once  the  % ^  are  chosen,  the  values 
of  <fi  and  all  the  sines  can  be  determined;  also,  functions  such  as  a.  X 

need  to  be  computed  only  NM  x  NK  x  NL  times.  Finally,  note  that  the  argument 
of  the  complex  Bessel  function  In  depends  only  on  the  two  indices  m  and  k  . 


L 


6d  J 
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Since  this  argument  is  independent  of  the  order,  ,  of  the  Bessel  functions, 
one  can  make  use  of  a  very  efficient  algorithm  based  on  the  recursion  relation 
(Ref.  261  to  calculate  all  orders  of  the  function  at  once,  for  a  given  m  and 
K  .  This  part  of  the  calculation  need  only  be  performed  NM  x  NK  times. 


Results  computed  from  Eqn.  (C-6)  were  compared  against  those  obtained 
using  the  numerical  integration  scheme  of  O'Hara  and  Smith“  for  several 
hundred  representative  integrals.  The  analytical  evaluation  proved  much 
faster,  and  agreement  within  the  1°>  accuracy  required  of  the  quadrature  calcu¬ 
lation  was  exhibited  in  the  great  majority  of  cases.  However,  for  the  higher 
values  of  m  and  *  ,  the  analytical  evaluation  gave  results  which  were  grossly 
in  error.  This  was  traced  to  the  fact  that  for  large  values  of  <1  and 
negative,  the  left  side  of  Eqn.  (C-4)  decays  exponentially  with  CL,  while 

■•j 

each  term  in  the  series  on  the  right  side  is  growing  exponentially.  This 
leads  to  similar  behavior  in  Eqn.  (C-61 ,  resulting  in  excessive  round-off 
error.  The  calculations  were  done  on  an  IBM  360/65  in  double  precision  arith¬ 
metic,  which  can  be  expected  to  yield  about  16  significant  digits.  For  these 
conditions  it  was  found  that  round-off  error  became  a  significant  problem  when 
CL  £  11,  in  which  case  the  program  automatically  reverts  to  the  numerical 
integration  method  cited  above.  Unfortunately,  when  a  is  large  b  is  also,  and 
the  resulting  highly  oscillatory  integrand  again  leads  to  a  time-consuming 
calculation.  Clearly  an  asymptotic  evaluation  of  the  integrals  in  this  limit 
would  be  preferable;  work  along  these  lines  has  been  initiated. 


For  l  i  2  ,  there  is  yet  a  further  means  of  streamlining  the  axial 
integral  evaluations.  This  is  based  on  the  following  easily  proven  symmetry 
propert  ies : 


C*C-xj  -  TZ/2  -  Czcx) 
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U'-~al 


C* C-  =  <-//*  C\X) 
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(C-~h) 


'3k 


C-i) 


.  t  *•  t  -l 

■I)  c 


8k 
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i  i  2  ,  m  c 


(C-"c 1 
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(-r )  5 


>,  2  ,  ?  > 


!  C  -  "’ll 
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*v-»  d  H 


C-  It) 


**•  (3< 


t*) 


cL<jMC*  » <5)  -  S'„ak(X-o)  -  <3  i  **«"  iC-'e 

It"  the  chordwi.se  collocation  points,  X(  .  are  symmetrical ly  disposed  about 
the  mid-chord,  X  *0,  then  use  ot"  the  above  relations  can  cut  computing  time 
and  storage  almost  in  half. 
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Figure  4  COMPARISON  OF  CHORDWISE  LOAD  DISTRIBUTIONS  AT  MID-ANNULUS 


WITH  ISOLATED  FLAT  PLATE  DISTRIBUTION 
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Figure  11  CONVERGENCE  OF  CALCULATIONS  VS.  NUMBER  OF  AZIMUTHAL 
MODES  FOR  CASES  1  AND  4 
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SOUND  LEVEL  CALIBRATOR 
GR  TYPE  1562- A 


Figure  15  ACOUSTIC  TEST  INSTRUMENTATION 
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b)  ROTOR  RPM  1710.5  OSCILLATOR  FREQUENCY  -  1312  Hi 

CALCULATED  BLADE  PASSAGE  FREQUENCY  1311  Hz 


Figure  16  OSCILLOSCOPE  RECORDS  OF  NOISE  SIGNAL  FROM  ROTOR  STATOR  INTERACTION 
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Figure  17  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND  PRESSURE 
LEVEL  FROM  ROTOR-STATOR  INTERACTION 
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Figure  17  (Cont.)  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND 
PRESSURE  LEVEL  FROM  ROTOR  STATOR  INTERACTION 
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Figure  17  (Cont.)  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND 
PRESSURE  LEVEL  FROM  ROTOR  STATOR  INTERACTION 
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Figure  18  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND  PRESSURE 
LEVEL  FROM  ROTOR  STATOR  INTERACTION 
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Figure  18  (Cont.)  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND 
PRESSURE  LEVEL  FROM  ROTOR  STATOR  INTERACTION 
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Figure  18  (Cont.)  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND 
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d)  ROTOR  RPM  -  1716 


Figure  18  (Cont.)  UNCORRECTED  1/3  OCTAVE  SPECTRUM  OF  OUTER  WALL  SOUND 
PRESSURE  LEVEL  FROM  ROTOR  STATOR  INTERACTION 


19  BACKGROUND  NOISE  CORRECTION  APPLIED  TO  MEASURED 
SOUND  PRESSURE  LEVELS 


MEASURED  FREQUENCY  RESPONSE  OF  MULTIFILTERS  IN  REAL  TIME  ANALYZER 


Figure  21  FREQUENCY  RESPONSE  OF  CONDENSER  MICROPHONE 


SOUND  PRESSURE  LEVEL  CORRECTED  FOR  BACKGROUND. 


Figure  23  RMS  WALL  PRESSURE  AT  SECOND  HARMONIC  OF  BLADE  PASSAGE 
FREQUENCY.  STATOR  STAGGER  ANGLE,  6  SM  =  37.2  deg 
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Figure  26  RMS  WALL  PRESSURE  AT  BLADE  PASSAGE  FREQUENCY.  STATOR 
STAGGER  ANGLE  A  SM  -  28.2  deg 
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SOUND  PRESSURE  LEVEL  CORRECTED  FOR  BACKGROUND, 
dB  re  2  x  10"*  dyrm/em^ 


Figure  27 


ROTOR  RPM 


RMS  WALL  PRESSURE  AT  SECOND  AND  THIRD  HARMONICS  OF  BLADE 
PASSAGE  FREQUENCY,  STATOR  STAGGER  ANGLE.  6  SM  -  28.2  deg 
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Figure  29  RMS  WALL  PRESSURE  AT  FIFTH  HARMONIC  OF  BLADE  PASSAGE 
FREQUENCY,  STATOR  STAGGER  ANGLE,  <5  =  28.2  deg 
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EXPERIMENT 


ROTOR  RPM 


•)  BLADE  PASSAGE  FREQUENCY 


Figure  30  COMPARISON  OF  THEORETICALLY  PREDICTED  SPL  WITH  EXPERIMENT 
STATOR  STAGGER  ANGLE  =  28.2° 


150 


{  uio/souAp  01  *  Z  J3M  0P)  3dS 


123 


STATOR  STAGGER  ANGLE 


EXPERIMENT 


ROTOR  RPM 


a)  BLADE  PASSAGE  FREQUENCY 


Figure  31  COMPARISON  OF  THEORETICALLY  PREDICTED  SPL  WITH  EXPERIMENT 
STATOR  STAGGER  ANGLE  =  3/.2° 
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